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This  report  contains  the  information  for  the  development  of  a  model  to  predict 
the  probability  of  explosion  propagation  between  adjoining  explosive  items, 
namely,  high  explosive  projectiles.  The  development  of  this  model  is  based 
on  the  sensitivity  of  an  encased  charge  to  primary  fragment  impciet. 

The  original  work  done  under  contract  to  Arthur  D.  Little,  Inc.,  was  extended 
to  provide  a  method  of  predicting  the  effectiveness  of  shielding  on  (cont) 
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projectile  spacing  requirements  and  also  to  predict  the  probability  of 
detonation  at  a  higher  confidence  level.  The  potential  practical  use  of  the 
model  is  judged  to  be  adequate  to  warrant  further  investigation.  When  used 
in  a  generalized  way,  the  model  can  be  a  valuable  tool  in  the  design  and 
layout  of  facilities  for  the  manufacture  and  storage  of  explosive  items. 
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SUMMARY 


This  report  describes  a  model  to  predict  the  probability  of 
explosion  propagation  between  adjoining  explosive  items.  The 
model  was  first  developed  by  Picatinny  Arsenal  under  contract  to 
Arthur  D.  Little,  Inc.,  of  Cambridge,  Massachusetts.  However, 
the  reliability  of  that  model  was  lower  than  the  desired  level. 
After  certain  modifications  had  been  made,  the  model  was  improved 
to  predict  probabilities  at  90  to  95  percent  confidence  level. 
The  model  was  further  extended  to  consider  cases  when  projectiles 
were  grouped  together,  funnels  attached  to  the  noses  of  the  pro¬ 
jectiles  and  shields  placed  between  them. 

The  development  of  the  model  was  centered  on  the  assumption 
that  the  sensitivity  of  a  projectile  to  detonation  depended  on 
the  following  properties  of  the  striking  fragment: 

1.  Velocity 

2.  Presented  area 

3.  Angle  of  impact 

4.  Amount  of  charge  in  projectile 

5.  Casing  tliickness  of  the  acceptor. 

The  model  was  finally  used  to  predict  the  probability  of 
detonation  propagation  between  individual  81mm,  105mm,  155mm  and 
8-in  projectiles,  and  also  between  groups  of  projectiles  on 
pallets.  Some  of  the  results  obtained  using  the  models, 
especially  for  those  cases  involving  individual  projectiles, 
compared  favorably  with  test  results.  In  the  case  of  grouped 
projectiles,  parameters  designed  for  single  and  groups  of  four 
projectiles  were  used  in  the  model  to  predict  the  probability  of 
detonation  propagation  between  two  similar  groups  of  sixteen  and 
seventy-two  105mm  and  81mm  projectiles,  respectively.  This 
approach  was  taken  merely  to  determine  the  adequacy  of  the  model. 
Further  tests  have  to  be  carried  out  to  determine  the  values  of 
the  parameters  in  question. 
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INTKOUUCTION 


Background 

An  analytical  model  to  predict  the  spacings  between 
adjoining  projectiles  necessary  to  achieve  any  given  level  of 
probability  that  a  detonation  will  not  propagate  from  one  to 
another  was  developed  by  Arthur  D.  Little,  Inc.  (ref.  1). 

This  model  is  based  on  theoretical  and  empirical  information 
that  has  been  compiled  by  previous  investigators  (refs.  2  to  6). 
A  limited  number  of  experiments  were  then  performed  by  ADL  to 
test  some  assumptions  made  and  to  evaluate  the  overall  adequacy 
of  the  model . 

The  development  of  the  model  was  centered  on  the  assumption 
that  the  sensitivity  of  an  acceptor  projectile  to  detonation  (50 
percent  of  the  time)  depends  on  the  following  parameters: 

1.  Presented  area  of  fragment 

2.  Velocity  of  fragment 

3.  Obliquity  angle  at  impact 

4.  Casing  thickness  of  projectile 

5.  Type  and  size  of  charge  in  projectile. 

Using  this  model,  safe  separation  distances  between 
projectiles,  which  will  result  in  the  detonation  being  propagated 
50  percent  of  the  time,  were  calculated  for  the  81mm,  105mm  and 
8-in  projectiles,  each  filled  with  Composition  B.  These  spacings 
were  compared  with  test  data  obtained  from  safe  separation  from 
tests  performed,  and  a  good  correlation  was  obtained  between  the 
two  sets  of  data  after  some  necessary  refinements  had  been  made. 

Objecti ves 

The  overall  purpose  of  this  study  is  to  develop  a  more 
reliable  model  to  predict  safe  separation  distances  between 
explosive  items. 

The  objectives  of  this  report  can  be  summarized  as  follows: 
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1.  To  review  the  ADL  model  and  incorporate  modifications 

in  it  so  as  to  extend  its  applicability  to  various 

practical  situations  common  to  AAP  facilities. 

2.  To  apply  the  modified  model  to  a  series  of  explosive 

propagation  cases  of  such  current  interest  as 

a.  Four  (2  x  2)  155mm  projectiles 

b.  Individual  8-inch  projectiles 

c.  Twenty-four  (6  x  4)  155mm  projectiles 

d.  Sixteen  (4  x  4)  105mm  projectiles 

and  to  compare  the  results  with  available  experimental  data. 

3.  To  increase  the  level  of  reliability  of  the  present 

(ADL)  model  and  extend  its  applicability  to  include  the 
following  situations: 

a.  Multiple  projectiles  -  The  probability  of 
explosive  propagation  between  pallets  loaded  with 
projecti 1 es . 

b.  Shielding  between  pallets. 

c.  Loading  funnels  -  Funnel  filled  with  Composition  B 

attached  to  the  nose  of  e?cn  projectile. 

4.  To  pinpoint  areas  where  additional  information, 

theoretical  and/or  experimental,  is  required  in  order 
to  further  refine  the  model. 

Format  and  Scope  of  Report 

The  entire  report  is  divided  into  two  parts  which  are 
further  subdivided  into  three  sections  each.  The  first  part  of 
the  report  describes  the  modified  (and  more  reliable)  model, 
considering  the  pertinent  topics  of  fragmentation  and 
sensitivity.  These  include  initial  fragment  velocity,  fragment 
mass  distribution,  effect  of  fragment  shape  on  velocity, 
sensitivity  of  encased  charge  and  the  effect  of  the  angle  of 
obliquity  on  the  sensitivity  of  the  charge.  The  probability  of 
detonation  for  a  given  distance,  a  procedure  for  calculating  this 
probability,  and  a  discussion  of  the  results  obtained  by  such  a 
procedure  are  also  presented  in  the  final  section  of  the  first 
part  of  the  report. 
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The  second  part  of  the  report,  which  begins  under  the 
section  "GROUPED  AND  SHIELDED  PROJECTILES",  deals  basically  with 
the  extension  of  the  model  (described  in  the  first  part  of  the 
report),  to  include  such  factors  as  the  number  and  velocity  of 
fragments  emitted  from  a  group-  of  projectiles.  The  first  few 
sub-sections  deal  with  the  description  and  the  development  of  the 
model  that  had  been  modified  to  predict  the  probability  of 
explosion  propagation  between  grouped  and/or  shielded 
projectiles.  The  final  section  is  a  brief  discussion  of  the 
results  that  can  be  obtained  from  the  new  model.  In  addition, 
appendix  A  contains  sample  calculations  incorporating  ideas 
developed  in  both  parts  of  the  study.  Appendix  B  provides  the 

format  of  the  input  data  cards  and  also  the  program  listing  and, 
finally,  appendix  C  presents  recommended  steps  to  follow  when 
projectiles  are  grouped  in  a  configuration  other  than  that  dealt 
with  in  this  study. 
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DESCRIPTION  OF  THE  MODEL 


General 

The  evaluation  of  the  adequacy  of  the  ADL  model  was  achieved 
through  tests  that  were  performed  by  James  Dobbie  et  al  (ref.  1). 
Such  tests  were  not  possible  for  this  study  and,  as  such,  all 
modifications  and  extensions  to  the  ADL  model  were  achieved 
through  the  analysis  of  experimental  and  theoretical  information 
collected  from  different  sources.  In  the  following  paragraphs, 
the  ADL  model  is  reviewed  and,  where  warranted,  changes  are  made 
accordingly. 

In  the  course  of  this  study,  it  was  necessary  to  develop  a 
computer  program  to  automate  the  prediction  calculation  at  each 
stage  of  its  development.  The  program  was  written  to  run  on  the 
IBM  1130  computer  system  and,  since  the  programming  of  the 
analytical  procedure  is  not  a  difficult  task,  the  program  can  be 
modified  to  run  on  any  system. 

Areas  of  Modification 

An  alternate  failure  criterion  as  expressed  by  the  Feist 
equation  (ref.  7)  was  defined  by  Picatinny  Arsenal.  The 
penetration  of  the  steel  casing  of  the  projectile,  even  without 
detonation,  was  a  mode  of  failure  that  had  to  be  considered.  As 
a  result,  the  sensitivity  equation  and  factors  affecting  the 
critical  velocity  had  to  be  re-examined. 

The  influence  of  certain  assumptions  and  coefficients 
(inherent  in  various  empirical  relationships  for  fragment 
characteristics  in  the  ADL  model)  on  the  predicted  separation 
distances  was  examined  and  alternate  formulations  were 
considered.  For  example,  the  value  for  Gurney's  constant  for 
Composition  B  used  in  the  ADL  model  is  approximately  11  percent 
higher  than  the  value  cited  in  reference  6. 

Fragment  Characteristics 

1.  Initial  Fragment  Velocity 

The  average  initial  velocity  of  fragments  from  an 
exploding  container  can  be  estimated  using  the  Gurney 
equation  which  assumes  that  the  charge  consists  of  an 
evenly  distributed  explosive  in  a  cylindrical  metal 
case : 

Vo  =  (2E)1/2[(C/M)/(1  +  C/2M)]1/2  (1) 
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where  (2E)l/2  =  8,800  fps  (Gurney's  constant) 

C  =  weight  of  charge 
M  =  weight  of  casing 

Vq  =  initial  velocity  at  center  polar  zone. 

The  C/M  ratio  takes  into  account  only  the  metal  in  the 
walls  of  the  casing  (ref.  6).  Therefore,  this  equation 
is  applicable  to  the  cylindrical  segment  of  the 
projectile  (the  sidewalls),  and  the  nose  and  base  area 
have  to  be  excluded.  Bearing  this  in  mind,  the  above 
equation  can  be  rewritten  as: 

Vo  =  (2E)1/2[2Q/(2  +  Q)]1/2  (2) 

where  Q  =  C/M  =  [pc(D  -  2x)2/[p^4(D-x)]  (3) 

D  =  outside  diameter  of  casing  (inches) 

pQ  and  Pfi,  =  densities  of  charge  and  metal  , 
respecti vely 

X  =  thickness  of  casing  (inches). 

Tests  performed  by  Gurney  with  TNT-filled  projectiles 
showed  a  fairly  good  correlation  with  the  formula  when 
(2E)1/2  =  8,000  ft/sec  was  used.  It  is  generally 
agreed  that  Composition  B  has  a  higher  density  than  TNT 
and  since  Gurney's  constant  varies  almost  linearly  with 
the  density  of  the  explosive,  some  researchers  have 
used  (2E)1/^2  =  8,800  fps  for  Composition  B.  As  will  be 
observed  later,  this  value  gives  a  better  result  than 
that  quoted  in  reference  6. 

2.  Variations  of  Fragment  Velocity  and  Mass  with  Polar 
Angl  e 

When  a  projectile  detonates,  fragments  are  projected  in 
many  directions  and  at  different  velocities. 
Experiments  have  shown  that  for  cylindrical 
projectiles,  the  greatest  density  of  fragments  is 
contained  in  a  narrov;  beamspray  generally  located  in 
the  central  polar  zone.  The  fragments  with  greatest 
velocities  are  also  found  in  this  region  which  is 
centered  at  the  90-degree  mark  from  the  nose  of  the 
projectile. 
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Tests  done  with  TNT-filled  155mm  projectiles  show  this 
narrow  region  to  be  centered  about  the  95-degree  mark, 
reinforcing  the  idea  that  the  beamspray  is  generally 
slightly  displaced  from  the  90-degree  mark  away  from 
the  point  of  initiation.  Using  these  test  results,  a 
formulation  for  the  initial  velocity  of  fragments  as  a 
function  of  the  polar  zone  was  achieved: 


Vi  =  VoFv(0) 

(4) 

where  Fv(0)  =  0.6474  -  0.026360  + 

0.000609502 

-  (3.08  X  10-6)03 

(5) 

50  j<0  £  950 

For  950  <0  <  1850,  (1900  _  0)  is 
equation  above. 

substituted  in 

the 

The  same  formulation  can  be  obtained  for  mass  of 
fragments.  Again,  using  the  same  test  results,  the 

weight  of  fragments  in  any  polar  zone  centered  at  0 
degrees  from  the  nose  of  the  projectile  can  be 
approximated  as: 

M(0)  =  CFm(0) 

(6) 

where  8^(0)  =  2.74326  -  0.0948750 

+  0.0009769402 

-  (2.388  X  10-6)03 

(7) 

500  <  0  <  950 

and  similarly,  950  <  0  _<  1850, 

equation  above. 

use  (1900  _  0) 

in 

C  =  total  weight  of  sidewalls  of  casing  (oz). 


To  determine  the  average  weight  of  fragments  ejected 
from  a  fragment  projectile,  the  equation  developed  by 
R.I.  Mott  can  be  used  (ref.  8): 

LnCN(m)]  =  LnCC'MA]  -  M/M^  (8) 

where  M  =  m^/^  (oz) 

C  =  C/2M/\3 
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M/\  =  is  a  fragment  distribution  parameter 
expressed  as 

Bx5/6(D  -  2x)1/3[1  +  x/(D  -  2x)]  (9) 

N(n))  =  no.  of  fragments  with  weight  >  rn 

B  =  explosive  constant  (approximately 
0.283  for  Composition  B) 

m  =  weight  of  fragment  (oz) 

X  =  casing  tiiickness  (in) 

D  =  average  outside  dianieter  of  casing  (in) 

C  =  casi ng  wei ght  (02) . 

It  should  be  noted  that  Mott's  equation  assumes  the 
projectile  casing  to  be  a  cylinder  with  a  uniform 
casing  thickness.  Although  this  is  hardly  the  case, 
satisfactory  results  have  been  obtained  with  this 
equation  (ref.  3). 

Mott's  equation  can  be  re-written  as: 

ln[N(m)]  =  ln[C/2M/\2]  - 

ln[2MA2N(m)/C]  = 

For  total  number  of  fragments,  m  =  0  and  therefore 
N  =  C/2Ma2 

Average  fragment  weight 

■  mo  =  C/N  =  2Ma2  (lo) 

The  ADL  expression  for  average  fragment  weight  is  based 
on  Heppner's  formulation  (ref.  4).  This  formulation 
(mg  =  2(KoDd/Vo)2)  gives  a  value  for  average  mass  of 
fragment  (for  tne  155mm  projectile)  that  is  40  percent 
higher  than  that  predicted  by  the  Mott  equation  which 
compares  well  with  test  results  (ref.  9). 
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3. 


Effect  of  Fragment  Shape  on  Velocity 


The  velocity  of  a  fragment  of  known  mass  and  presented 
area  at  a  distance  R  feet  from  point  of  initiation  can 
be  approximated  by: 


Vs 

where  k' 
Vs 

Vo 

R 

A 

Pa 

Cd 


Voe(-k'R/m^/^)  (ref.  10)  (11) 

(A/m2/3)p^C[) 

striking  velocity  of  fragment  at  distance 
R  (ft) 

initial  velocity  of  fragment  (fps) 

distance  from  source  (ft) 

presented  area  of  fragment  (in^) 
or  fragment  impact  area 

air  density  (oz/in3) 

drag  coefficient. 


The  above  equation  can  be  rewritten  as: 

Vs  “  Voe("0*0025Sm  ^Z^)  (ref.  1)  (12) 

S  =  distance  from  donor  source  (in) 


m  =  mass  of  fragment  (grain) 


Sensiti vity 

The  sensitivity  equation  given  by  ADL  is  based  on  the 
original  work  of  Slade  and  Dewey  (ref.  11)  in  which  cylindrical 
fragments  were  fired  against  bare  and  cased  charges  (tetryl  and 
Composition  B).  It  was  concluded  that  the  velocity  for  50 
percent  initiation  was  a  function  of  contact  area  and  not  the 
mass  of  the  striking  fragment,  and  this  velocity  increased 
approximately  linearly  with  the  thickness  of  the  cover  plate  of 
charge.  Experiments  carried  out  by  Picatinny  Arsenal  (ref.  12) 
tend  to  disprove  Slade  and  Dewey's  theory.  However,  enough  data 
was  not  available  to  reach  a  decisive  conclusion. 

To  be  able  to  design  the  model  to  predict  the  50  percent 
velocity,  ADL  used  the  same  criteria  used  by  Slade  and  Dewey; 
that  is,  the  50  percent  velocity  which  was  defined  as  the 
velocity  at  which  break-up  and  total  detonation  of  charge  were 
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equally  probable.  This  velocity  was  calculated  whenever  possible 
by  taking  the  arithmetic  mean  of  the  highest  velocity  for 
break-up  and  the  lowest  for  detonation. 

Since  it  is  an  objective  of  this  study  to  improve  the 
reliability  of  the  model,  it  was  necessary  to  re-analyze  the 
Slade  and  Dewey  data  using  the  lowest  velocities  at  which  some 
deflagration  occurred.  By  considering  these  velocities  as 
critical,  and  also  considering  the  linear  variation  of  these 
velocities  with  the  thicknesses  of  the  cover  plate,  a  formulation 
was  achieved  defining  the  velocity  required  for  partial 
detonation  (see  figs.  1  and  2). 

For  a  cylindrical  fragment  of  diameter  d  striking  at  zero 
impact  angle,  the  required  velocity  can  be  expressed  as; 

Vs  =  (1  +  Kix/d)Jd-3  (13) 

(1820.74  Composition  B 
J  =  ( 

(1030.77  tetryl 

(0.61  Composition  B 
a  =  ( 

(0.42  tetryl 

Vs  =  striking  velocity  (fps) 

X  =  thickness  of  acceptor  plate  (in) 
d  =  diameter  of  fragment  (in) 

=  ADL  constant  (0.6  -  0.8) 

J  and  a  are  proportionality  constants 

Velocity  predicted  by  the  above  formula  is  less  than  that 
predicted  by  ADL's  sensitivity  equation,  but  greater  than  the 
velocity  required  to  just  penetrate  a  steel  cover  plate  of 
thickness  x  (in). 

Effect  of  Obliquity 

Analyses  of  data  from  Slade  and  Dewey  showed  that  the  lowest 
velocity  required  for  partial  detonation  varied  approximately  the 
same  way  as  the  50  percent  velocity  varies  with  the  angle  of 
impact  of  the  fragment.  Therefore,  any  analysis  here  would  be 
identical  to  that  presented  by  ADL  (see  fig.  3). 
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FRAGMENT  VELOCITY  (meters  per  second) 


FIGURE  3.  VARIATION  OF  VELOCITY 
WITH  IMPACT  ANGLE 

(from  SloSt  ft  Dt«ty  ,  BRL  Rtporf  No.  1021  ) 
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COMPUTATION  OF  PROBABILITY  OF  DETONATION  FOR  A  GIVEN  SPACING 


Probability  of  Detonation  as  a  Function  of  Distance 

From  information  available  in  ADL's  report,  we  can  estimate 
the  expected  number  of  fragments  in  any  zone.  We  can  also 
calculate  the  expected  number  of  effective  hits  (i.e.,  hits 
capable  of  causing  detonation)  at  an  acceptor  at  some  distance  S 


from  expl osi on : 

Ni  =  N-j  f -j  Fqci  F^i 

(14) 

wliere  F^-j  =  (D-j  -  3x-j)/2  S 

(15) 

fi  =  k3expC-(mi)^/2/M;\] 

(16) 

”  CF|'n(0)/mo 

(17) 

D-j  =  average  outside  diameter  of  acceptor  (in) 


x-j  =  average  case  thickness  of  acceptor  (in) 

m-j  =  critical  mass  in  polar  zone 

Fzi  =  fraction  of  the  i^h  lO-degree  zone 

Ma  =  fragment  distribution  parameter  (eq.  9) 

Fccj  =  fraction  of  the  360-degree  azimuthal 
angle  intercepted  by  the  acceptor 

f-j  =  fraction  of  fragments  that  strike 
acceptor  that  are  supercritical 

N-j  =  expected  number  of  fragments  in  the 
i^^  zone 

N-j  =  expected  number  of  supercritical  hits 
in  the  i^^  zone 

The  total  number  of  supercritical  hits  is 
zone  _ 

N  =  ^  Ni  (18) 

i=l 
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and  the  chance  of  at  least  one  supercritical  hit  at  the  acceptor 
is 


P  =  1  -  (19) 

The  relationships  discussed  above  supply  us  with  a  method  for 
predicting  the  probability  of  detonation  propagation  for  a  given 
distance.  The  derivation  of  these  equations  was  done  by  ADL. 
There  was  no  need  to  repeat  the  derivations  here,  and  since  this 
study  is  an  extension  of  that  done  by  ADL,  the  interested  party 
should  consult  the  ADL  report. 

Procedure  for  Computing  the  Probability  of  Detonation  as  a 
Function  of  Spacing 

All  the  equations  presented  in  the  earlier  sections  have 
been  collected  here  in  order  to  design  a  procedure  for 
calculating  the  probability  of  detonation.  The  derivation  of 
some  of  these  equations  can  be  found  in  the  report  prepared  by 
Arthur  D.  Little,  Inc.  The  probability  of  detonation  of  acceptor 
shall  be  calculated  from  the  following  steps: 

Step  1: 

Measure  from  a  scaled  drawing  of  the  projectile: 

D  =  the  average  outside  diameter  at  the  middle  section 
X  =  the  average  casing  thickness. 

Step  2: 

Calculate  the  initial  velocity  of  fragments  using  equation 

(2): 


Vo  =  (2E)1/2[2Q/(2  +  Q)]1/2. 

For  this  velocity,  the  charge-to-metal  ratio  has  to  be 
calculated.  This  is  achieved  by  using  equation  (3). 
Gurney's  constant  should  be  taken  as  8,880  fps  for 
Composition  B. 

Step  3: 

On  a  scaled  drawing  as  shown  in  figure  A.l,  with  donor  and 
acceptor  projectiles  spaced  S  (in)  center  to  center  and 
their  bases  on  the  same  level,  lay  off  the  dividing  rays  for 
polar  zones  of  10  degrees  centered  at  integral  multiples  of 
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10  degrees.  Determine  the  zones  in  which  fragments  can  hit 
tlie  acceptor,  ignoring  the  base  plate  and  fuze  section. 
Estimate  the  fraction  F2i  of  the  zone  over  wliich  the 
vulnerable  part  of  the  acceptor  extends.  F^-j  =  1.0  for 
zones  completely  spanned.  If  S  is  large  enough  that  only 
the  central  zone  is  involved,  then 

Fz  =  ISH/ttS 

wliere  H  =  height  of  vulnerable  portion  of  projectile.  Also, 
estimate  the  average  outside  diameter  (in)  and  thickness  of 
casing  (in)  in  each  zone. 

Step  4: 

For  each  zone,  compute  the  critical  mass  m-j  (grain)  by 
iteration  from  the  following  equations: 


Fp.  =  e(-0-0025Smi-l/3)  (20) 

Vi  =  VoFviFoi  (21) 

d  =  (4FAA)l/2(mi/k)l/3*  (^2) 

Vi  =  (1  +  KiX/d)Jd-a  (13) 


Starting  with  m-j  =  mQ,  calculate  Fp-j  and  V-j  using  equation 
(5)  to  calculate  F\/-j.  For  this  same  value  of  m-j,  calculate 
d  from  equation  (22)  above,  and  finally  calculate  V-j  using 
equation  (13).  Repeat  this  until  the  V-j's  given  by 
equations  (21)  and  (13)  do  not  differ  significantly.  F/\, 
tlie  fraction  of  the  contact  area  of  the  fragment  that  is 
effective,  is  taken  as  0.50. 

Step  5: 

For  each  zone,  compute  the  total  number  N-j  of  fragments,  the 
azimuthal  factor  F  -j ,  and  the  fraction  f-j  of  the  striking 
fragments  that  have  mass  exceeding  the  critical  mass  m-j 
using  the  following  operations: 


*Equation  (22)  was  derived  by  ADL  for  the  sensitivity  equation  to 
take  into  account  the  mass  of  the  striking  fragment.  The 
iterative  process  in  this  section  can  be  time-consuming  unless 
one  of  the  advanced  numerical  techniques,  for  example,  the  secant 
method, is  used. 
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Ni  =  CF^(0)/mo 

(6) 

=  (Di  -  2Xi)/2TTS 

(15) 

fi  =  Fie-(2mi/mo)^/^ 

F(y,(0)  is  calculated  from  equation  (7).  FI,  another  factor 
derived  by  ADL,  considers  the  effect  of  non-normal  impacts 
and  it  is  estimated  to  be  0.2  for  Composition  B.  C  is  the 
total  weight  of  the  casing  sidewalls  (grain). 

Step  6: 

The  probability  of  detonation  can  now  be  calculated  using 
the  following: 


z 


N  =  E  NjF  jFzifi 

(18) 

II 

1— * 

P  =  1  -  e(-f^) 

(19) 

z  =  number  of  zones  involved. 

Discussion  of  Results 

To  improve  the  reliability  of  the  ADL  model,  calculations 
were  done  for  the  81mm  (M374A1  Comp  B)  projectiles  and  results 
compared  with  those  presented  in  re‘’'erence  13.  The  calculations 
were  performed  on  the  computer  and  the  input  and  output  data  are 
presented  in  appendix  A. 

Certain  notable  changes  made  in  the  ADL  model  were: 

1.  The  value  of  the  shape  factor  k  (or  ballistic  density): 
This  value  expresses  the  relationship  between  fragment 
mass  and  projected  areas.  A  k  factor  of  660  grains  per 
cubic  inch  was  found  to  yield  a  more  accurate  value  of 
the  separation  distance  than  that  quoted  in  the  ADL 
report  (640  grams/cu  in).  Reference  9  makes  use  of  tlie 
same  value. 

2.  Area  Factor  F/\:  Calculations  sliowed  that  a  value  of 
0.5  provided  better  results  than  0.75  as  presented  in 
the  ADL  report.  Tests  stiowed  the  average  fragments  had 
at  least  one  pointed  side,  thus  reducing  the  ratio  of 
the  effective  area  to  the  presented  area. 
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3.  Protection  Coefficient  K^;  A  closer  examination  of  the 
results  presented  in  reference  3  indicated  =  0.8  to 
be  a  better  approximation. 

With  the  above  changes  made,  the  necessary  parameters  and 
dimensions  of  the  projectile  were  inputed  into  the  program  in  the 
format  shown  in  appendix  A.  For  a  distance  of  21.08  inches 

between  individual  81mm  projectiles,  the  probability  of 
detonation  propagating  from  one  to  the  other  was  found  to  be  6.4 
percent.  Tests  performed  by  Picatinny  Arsenal  (ref.  13)  gave  a 

probability  of  5.1  at  95  percent  confidence  level  for  the  same 
separation  distance.  The  insignificant  difference  between  the 
values  predicted  by  the  updated  model  and  experiments  is  a 
consequence  of  the  conservative  assumptions  made  in  the  model. 

Calculations  were  also  performed  for  the  155mm  projectiles 
with  separation  distances  of  1.68  m  (66.1  in)  and  2.29  m  (90.1 

in).  For  both  distances,  the  model  predicted  a  probability  of 

detonation  of  12.7  and  8.7  percent,  respectively.  These 
predictions  were  compared  with  the  test  results  documented  in 
reference  14.  For  a  separation  distance  of  2.29  m  (90.1  in),  one 
out  of  12  projectiles  resulted  in  a  high  order  detonation  whicfi 
is  equivalent  to  a  probability  of  8.3  percent.  However,  for  a 
confidence  level  of  95  percent  and  11  observations  with  no 
reactions,  the  probability  of  explosion  propagation  will  fall 
between  zero  and  28.5  percent. 

For  a  separation  distance  of  1.68  m  (66.1  in)  and  a 
confidence  level  of  95  percent,  the  true  value  of  the  probability 
of  explosion  propagation  was  between  the  zero  to  45.9  percent 
interval.  This  relatively  wide  interval  is  based  on  six 
observations,  as  the  number  of  observations  increases,  the 
confidence  level  increases  (the  interval  decreases). 

Calculations  were  again  performed  for  the  same  projectiles 
and  separation  distance,  but  this  time  the  new  alternate  failure 
cri ter ioni(perforation  of  casing  even  without  high  order  detonation) 
was  considered.  The  model  predicted  a  probability  of  detonation 
of  36  percent  for  a  separation  distance  of  21.08  inches  between 
individual  81mm  projectiles.  The  accuracy  of  this  prediction 
could  not  be  determined  due  to  the  unavailability  of  test  data  in 
which  perforation  of  casing  was  considered  a  failure  criterion. 
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GROUPED  AND  SHIELDED  PROJECTILES 


General 

The  reliability  of  the  model  developed  by  ADL  was  improved 
in  part  1  of  this  report.  The  model  designed  to  predict  the 
probability  of  detonation  propagation  between  adjacent  explosive 
items  was  refined  to  incorporate  the  alternate  failure  criterion 
proposed  by  ARRADCOM. 

As  part  of  the  overall  objective  of  this  study,  the  model  is 
extended  to  predict  "safe  distances"  for  a  much  broader  range  of 
practical  configurations  such  as: 

1.  Four  (2  X  2)  155mm  projectiles  on  pallets 

2.  Twenty-four  (6  x  4)  155mm  projectiles  on  pallets 

3.  Sixteen  (4  x  4)  155mm  projectiles  on  pallets. 

Certain  elements  that  influence  the  safe  separation  distance 
are  examined.  These  include  situations  when  plate  or  bar  shields 
are  placed  between  stacked  projectiles  and  when  funnels  are 
attached  to  the  noses  of  the  projectiles  on  the  pallets. 

The  depth  and  accuracy  of  this  portion  of  the  study  is 
adversely  limited  by  the  unavailability  of  experimental  data. 
Most  of  the  equations  and  discussions  presented  in  this  section 
are  based  entirely  on  the  analysis  of  data  available  in  Technical 
Report  No.  3664,  "Fragment  Hazard  Program",  by  Richard  T.  Ramsey, 
et  al  (ref.  9). 

Characteristics  of  Fragments  Emitted  from  a  Group  of  Projectiles 

In  tlie  case  of  an  individual  projectile,  it  was  assumed  that 
the  source  of  tfie  emitted  fragments  was  rotational ly  symmetric 
about  the  longitudinal  axis  and  thus  the  properties  of  the 
fragments  were  functions  of  the  polar  angle  only.  Test  results 
show  that  this  assumption  is  not  valid  for  grouped  or  stacked 
projectiles.  Properties  of  fragments  emitted  from  such  a  group 
are  functions  of  the  polar  and  azimuthal  angles. 

However,  when  the  projectiles  are  grouped  in  a  "square 
matrix"  arrangement,  it  can  be  conservatively  assumed  that 
rotational  symmetry  for  such  a  group  exists.  Tests  documented  in 
reference  9  show  the  above  assumption' to  be  fairly  accurate. 
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Due  to  the  limited  availability  of  data  on  group 
projectiles,  consideration  would  be  given  to  the  square  matrix 
arrangement.  Extension  of  the  model  to  include  any  stack 
configuration  should  not  be  a  difficult  task. 


© 


Figure  5 


Other  assumptions  that  have  to  be  made  for  group  projectiles  in  a 
square  matrix  arrangement  include: 

1.  The  metal  in  the  shaded  region  in  figure  4  does  not 

enter  the  fragment  mass  distribution. 

2.  Targets  1,  2,  3  and  4,  situated  at  the  same  distance 

from  the  donor,  have  equal  probabilities  of  being  hit 
by  fragments  emitted  from  the  group  (this  follows  from 
the  assumption  that  the  source  is  rotational ly 

symmetric  about  its  longitudinal  axis). 

It  should  be  emphasized  that  the  above  assumptions  would 

lead  to  results  that  are  very  conservative  and  they  are  also  only 

valid  for  the  arrangement  shown  in  figure  5.  When  a  number  of 
projectiles  other  than  four  (2  x  2)  are  grouped  together  in  a 
configuration  other  than  that  described  above,  the  steps 
described  in  appendix  C  have  to  be  followed.  Again,  it  should  be 
emphasized  that  the  accuracy  of  these  steps  could  not  be 
determined  due  to  insufficient  experimental  data. 

Initial  Fragment  Velocity 

The  velocities  of  fragments  emitted  from  a  group  of  four  (2 
X  2)  projectiles  were  found  to  be  twice  those  from  an  individual 
projectile  between  the  60-  and  120-degree  polar  angles.  Outside 
this  zone,  the  fragment  velocities  were  approximately  1-1/2  times 
those  for  a  single  projectile  (ref.  9). 
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The  beamspray  was  also  found  to  be  centered  around  the 
103-degree  mark  unlike  the  95-degree  mark  in  the  case  of  a  single 
projectile.  This  8-degree  shift  affects  the  velocity  parameter 
Fy(0).  To  determine  Fy(0)  for  a  group  of  four  projectiles,  (0  - 
8)  should  be  used  in  equation  (5)  for  13  _<  0  _<  103  and  for  polar 
angles  between  103°  and  198°  (198  -  0)  should  be  used. 

Fragment  Mass  Distribution 

From  assumption  2,  the  distribution  parameter  F^(0)  defined 
in  equation  (7)  can  also  be  used.  However,  due  to  the  8-degree 
shift  in  the  location  of  the  beamspray,  (0  -  8)  should  be 
substituted  for  0  in  equation  (7)  for  13  _<  0  _<  103  and  (198  -  0) 
is  used  for  angles  outside  this  limit.  The  mass  of  fragments 
emitted  from  a  stack  of  projectiles  in  a  polar  zone  centered  at 
angle  0  can  be  expressed  as: 

M(0)  =  CAFni(0)  (6) 

where  Ca  is  the  fraction  of  the  total  weight  of  metal  involved  in 
the  mass  distribution. 

Figure  6  below  shows  how  this  value  can  be  estimated.  Test 
results  show  this  method  to  provide  estimates  that  are  within  the 
tolerance  limit  of  this  study. 


°  (azimuth) 


f80^ 


Figure  6. 


Tests  performed  with  155mm  projectiles  in  reference  9  showed 
also  that  the  total  weight  of  fragments  from  an  individual 
projectile,  collected  in  an  area  bounded  by  45-  to  135-degree 
polar  zone  was  61,378  grains.  The  average  weight  of  metal 
collected  in  a  similar  area  for  a  group  of  4  projectiles  was 
175,747  grains  (giving  a  ratio  of  approximately  3).  From  figure 
6,  the  model  shows  3/4  sectors  of  projectile  metal  available  in  a 
90°  azimuthal  sector.  This  gives  a  ratio  of  3  (3/4  projectile 
from  a  group  of  4  and  1/4  from  a  single  projectile). 
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The  slopes  of  the  curves  in  figure  7  show  that  the  average 
weight  of  the  fragments  emitted  from  a  single  and  group 
projectiles  was  the  same.  This  is  the  expected  result  as 
predicted  by  Mott's  equation  presented  in  the  first  part  of  this 
report  (eq.  10). 

Effect  of  Siiields 

1.  Plate  Shield 


Perforation  data  for  steel  impacting  on  several 
metallic  materials  including  steel  and  aluminum  was 
collected  and  analyzed  by  Ballistic  Researcit 
Laboratories  as  part  of  project  THOR  (ref.  2).  These 
test  results  showed  that  the  residual  velocity  of  a 
fragment  after  perforation  of  a  medium  could  be 
expressed  as; 


Vr  -  V3 

where 

^s 

t 

A 

ms 

e 


10C(tA)“mf (sec  e)^Vs  (23) 

fragment  residual  velocity  (fps) 
fragment  striking  velocity  (fps) 
thickness  of  target  (in) 
average  impact  area  (in^) 
initial  weight  of  fragment  (grain) 
obliquity  angle 


c,  oi,  B,  Y  and  X  are  constants  and  these  values  are 
listed  below  in  table  1  for  mild  homogeneous  steel  and 
Aluminum  Alloy  2024T-3. 


Table  1 


Target  Material 

c 

a 

Y 

A 

Mild  homogeneous  steel 

6.399 

0.889 

-0.945 

1.262 

0.019 

Hard  homogeneous  steel 

6.475 

0.889 

-0.945 

1.262 

0.019 

Aluminum  Alloy  2024T-3 

7.047 

1.029 

-0.072 

1.251 

-0.139 

Cast  iron 

4.840 

1.042 

-1.051 

1.028 

0.523 
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N(m),  NO.  OF  FRAGMENTS  WITH  MASS  GREATER  THAN 


lO.OOOr 


FIGURE  7  MOTT  PLOTS  FOR  ONE(I), 
FOUR(2x2)  a  EI6HT(4x2) 
GROUPS  OF  PROJECTILES 

(from  R.T  Romooy,  Frofmont  Hozord  Invootiflofion  Pro«rOfii  ) 
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By  setting  the  residual  velocity  V^-  =  U  in  equation 
(23),  a  good  approximation  of  the  protection  velocity 
can  be  obtained. 

Vq  =  lu''l(eA)“lm5^1(sec.  6  (24) 

v^here  Vq  =  protection  velocity  (fps) 

e,  A,  mg  and  6  have  been  defined  previously. 

c^,  a]^,  and  yi  are  constants  whose  values  are  listed 
i n  Table  2. 

Tfie  protection  velocity  Vq  is  defined  to  be  the  highest 
striking  velocity  below  the  ballistic  limit  for  which 
probability  of  perforating  a  given  target  is  zero. 


Table  2 


Target  Material 

Cl 

«1 

ei 

Yl 

Aluminum  Alloy  2024T-3 

6.185 

0.903 

-0.941 

1.098 

Mild  homogeneous  steel 

6.523 

0.906 

-0.963 

1.286 

Hard  homogeneous  steel 

6.601 

0.906 

-0.963 

1.286 

Cast  iron 

10.153 

2.186 

-2.204 

2.156 

After  perforation,  the  mass  of  tfie  fragment  varies  and 
again  an  estimate  of  the  residual  mass  of  the  fragment 
can  be  obtained  by: 

m,-  =  nig  -  10*^ (tA)“m|(sec  (25) 

Table  3  below  lists  the  value  of  the  constants  in 
equation  (25). 
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Table  3 


Target  Material_ c_ a_ g_ x ^ 

Mild  homogeneous  steel  -2.507 


Hard  homogeneous  steel  -2.264 
Aluminum  Alloy  2024T-3  -6.663 
Cast  iron  -9.703 


0.138 

0.835 

0.143 

0.661 

0.346 

0.629 

0.327 

0.880 

0.227 

0.694 

-0.361 

1.901 

0.162 

0.673 

2.091 

2.710 

The  critical  mass  (weight  of  fragment  to  cause 
detonation)  is  calculated  as  a  function  of  the  residual 
velocity  and  mass.  A  decrease  in  the  velocity  and 
weight  of  fragment  after  perforation  would  reduce  the 
probability  of  such  a  fragment  being  critical. 

2.  Bar  Shields 

Available  test  results  with  bar  (or  pipe)  shields  are 
not  sufficient  enough  to  be  able  to  establish  empirical 
relationships.  However,  tests  performed  by  C.  Anderson 
and  R.  Rindner  of  Picatinny  Arsenal  (ref.  14)  with 
steel  and  aluminum  rods  placed  between  single  and 
multiple  155mm  projectiles  (M107  Comp  B)  provided  the 
following  observations: 

a.  The  fragments  from  the  donor  projectile  which  are 
"head-on"  to  the  critical  region  (a  portion  of  the 
area  of  the  projectile  containing  the  charge)  of 
the  acceptor  are  deflected  by  the  bar  shield, 
while  those  fragments  which  miss  the  shield  would 
only  have  a  grazing  impact  on  the  projectile. 
However,  in  the  case  of  multiple  projectiles, 
these  grazing  fragments  might  strike  normal  to 
other  projectiles  in  the  group. 

b.  Aluminum,  besides  being  lighter  and  less 
expensive,  absorbs  more  energy  per  unit  weight 
during  deformation  (plastic)  and  behaved  more 
uni formly  than  steel . 

The  existence  of  shields  (bar  or  plate)  between  donor 
and  acceptor  projectiles  reduces  the  probability  of  an 
acceptor  projectile  being  hit  by  a  fragment.  A  very 
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Table  4.  Parameter  values 


Shape  Factor  (Ballistic  Density)  K  (grain/in^)  =  660.0 

Sensitivity  Parameter  B  =  2,140.0 

Area  Factor  FA  =  0.5 

Constant  J  (sensitivity)  =  1,820.74 

Constant  A  (sensitivity)  =  0.61 

Constant  B  (Mott's  equation)  =  0.283 

Average  Mass  Coefficient  KO  =  4,500.0 

Protection  Coefficient  K1  =0.8 

Impact  Angle  Factor  FI  =  0.2 

Gurney's  Constant  VC  (ft/sec)  =  8,800.0  for  Composition  B 

=  8,000.0  for  TNT 

Charge  Density  (Ib/in^)  =  0.058  for  Composition  B 

=  0.056  for  TNT 


Table  5.  Various  explosive  constants  used  in  calculations 


Explosi ve 
type 

Gurney' s 
Constant 
(fps) 

Expl osi ve 
Constant  B 
(eq.  9) 

Constant  J 
(eq.  13) 

Constant  A 
(eq.  3) 

TNT 

8,000.00 

0.300 

- 

- 

Composition  B 

8,800.00 

0.283 

1,820.74 

0.61 

Tetryl 

7,460.00 

0.240 

1,030.77 

0.42 
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conservative  estimate  of  such  a  probability  can  be 
obtained  from  the  model  developed  in  tliis  study. 

Effect  of  Funnels 

The  extra  v/eight  of  funnels  attached  to  the  noses  of 
projectiles  enter  the  fragment  weight  distribution.  The 
assumption  inherent  in  Gurney's  equation,  namely,  the  projectile 
casing  is  cylindrical,  can  no  longer  be  considered  valid  here. 
However,  if  the  following  assumptions  are  made,  the  effects  of 
funnels  can  be  dealt  with  but  the  validity  of  these  assumptions 
cannot  be  justified  due  to  insufficient  experimental  data. 

1.  The  extra  weights  of  the  casing  of  the  funnel  and 
enclosed  charge  are  taken  into  account  when  the  C/M 
ratio  (eq.  3)  is  calculated,  and  initial  velocity  is  a 
function  of  this  ratio  and  not  of  the  outside  diameter 
and  tliickness  of  casing. 

2.  Perforation  of  any  funnel  in  the  group  will  cause 
detonation  of  the  entire  group  of  projectiles. 
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DEVELOPMENT  OF  MODEL  TO  PREDICT  PROBABILITY 
OF  DETONATION  FOR  STACKED  PROJECTILES 


Probability  of  Detonation  as  a  Function  of  Distance  and  Shielding 

All  of  the  extensions  to  the  old  model  discussed  in  the 
previous  section  are  grouped  here  to  facilitate  the  development 
of  the  nev/  model.  The  probability  of  detonation  can  be  computed 
using  the  following  equations: 

1.  Fragment  initial  velocity: 

Vo  =  (2E)1/2[20(2  +  Q)]l/2.  (2) 

For  a  group  of  four  (2  x  2)  155mm  projectiles,  initial 
velocity  of  fragments  in  a  polar  zone  centered  at  0  is: 

V(0)iN  =  nVoFV(0)  (4) 


(2.0 

60  ^  0  ^  120 

where 

n  =  ( 

(1.5 

otherwi se 

where 

Fv(0) 

=  0.6474  -  0.02636(0  -  8) 

+  0.0006095(0  -  8)2 
-  (3.08  X  lO-6)(0  -  8)3  (5) 

for  130  <  0  <  1030. 

For  1030  ^  0  ^  198°,  use  (198  -  0)  in  the  equation 
above. 

2.  Number  of  fragments: 

ln[N(m)]  =  ln[C*MA]  -  M/M^  (8) 


where  M  =  ml/2  (oz) 

C*  =  Ca/2Ma3 

Ma  =  Bx5/6(d  -  2x)1/3[1  +  x/(D  -  2x)]  (9) 

The  notations  used  above  have  been  defined  previously. 
The  new  constant  Ca  was  also  defined  earlier. 
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3.  Average  mass  of  fragment: 


Total  number  of  fragments  involved  in  fragment  mass 
di stribution: 

N  =  Ca/2Ma2 

Average  fragment  weight:  mo  =  C/\/N  = 

4.  Striking  velocity  at  a  given  distance  from  donor: 


Vs  =  Voe(-K'R/ml/^) 

(11) 

=  Voe(-0-0025Sm-l/3) 

(12) 

where  S  =  distance  (in) 

m  =  mass  of  fragment. 

When  a  shield  is  present  between  donor  and  acceptor 
projectiles,  striking  velocities  at  surface  of  shield 
and  at  acceptor  should  be  calculated. 


For  a  plate  (bar)  shield,  striking  velocity  is: 


where 


Vsp  =  V,Ne(-“-0025Siiiisl/3)  (12) 

Vgp  =  striking  velocity  of  fragment  at  surface 
of  plate  shield  (■'"ps) 


=  initial  velocity  of  fragment  at  donor 
source  (fps) 


Si  =  distance  between  donor  and  surface  of 
shield  next  to  donor  (in) 


ms  =  initial  fragment  mass  (grain) 


Assuming  the  shield  is  perforated  by  the  fragment, 
striking  velocity  of  fragment  at  acceptor  surface 
i  s : 

Vsr  =  Vre(-0-0025S2nir“^/^)  (12) 

where  1/5^  =  striking  velocity  of  fragment  at 
acceptor  (fps) 
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Vp  =  residual  velocity  of  fragment  at 
sliield  (fps) 

$2  =  distance  between  shield  and 
acceptor  projectiles  (in) 

=  residual  weight  of  fragment  (oz) 

It  should  be  noted  that  Si  is  the  maximum  distance 
(in)  between  the  donor  pile  and  the  shield  that 
would  result  in  the  most  effective  use  of  the 
shield.  Figure  8  shows  how  this  distance  can  be 
calculated.  A  basic  assumption  here  is  that  any 
fragment  that  does  not  perforate  the  shield  is 
either  blocked  by  it  or  deflected  away  from  the 
acceptor  projectile. 


Donor  Shield  Acceptor 


Figure  8. 


5.  Minimum  striking  velocity  of  a  fragment  of  mass  (m)  to 
cause  detonation: 

Vs  =  (1  +  Kix/d)Jd-a  (13) 

d  =  (4FA/?r)l/2(mi/k)l/3  (22) 

v^here  F/\  =  constant  of  proportionality  between 

impact  and  effect  areas  of  the  fragment. 

The  striking  velocity  at  the  surface  of  the  acceptor 
projectiles  V51  is  calculated  as  a  function  of  the 
residual  velocity  after  perforation  of  shield,  taking 
into  account  the  drag  effects. 
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After  perforation  of  the  shield,  the  mass  of  the 
fragment  changes  and  so  does  its  diameter.  Equation 
(19)  should  be  changed  to: 

d  =  (4FATr)l/2(mir/K)l/3  (22) 

where  m-j,-  =  residual  mass  in  the  i'*^^  zone  (grain). 

The  equations  presented  here  would  now  be  used  to 
determine  the  probability  of  detonation  for  grouped 
projectiles. 

Procedure  for  Calculating  Probability  of  Detonation  Propagation 
Between  Stacked  Projectiles 

The  probability  of  detonation  propagation  between  stacked 
projectiles  (4  in  a  group)  can  be  determined  from  the  following 
steps : 

Step  1. 

From  a  scaled  drawing  of  a  single  projectile  (as  shown  in 
figure  A.l),  determine  the  following: 

D  =  outside  diameter  of  projectile  (in) 

X  =  average  casing  thickness  (in)  should  be  measured 
around  the  center  of  the  projectile. 

If  funnels  are  involved  in  the  stack,  then  determine: 

M  =  total  weight  of  casing  (oz).  This  includes  weight 
of  the  funnel  casing  and  the  sidewalls  of  the 
projectile. 

C  =  the  total  charge  weight  (oz)  both  in  projectile 
and  funnel . 

Step  2. 

Calculate  initial  velocity  of  fragments  from  a  single 
projectile  (with  or  without  funnel)  from  the  following 
equation : 

Vo  =  (2E)1/2[2Q/(2  +  Q)]1/2  (2) 
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where  Q  =  Pc{D  -  2x)2/[p^{D  -  x)] 

for  a  cylindrical  projectile  without  a  funnel 

or  where  Q  =  C/M 

for  a  projectile  with  a  funnel. 

Step  3.  (If  shields  have  to  be  used.) 

On  a  scaled  drawing  with  donor  and  acceptor  stacks  in  a 
horizontal  position,  space  the  distance  S,  and  determine  the 
maximum  distance  between  the  shield  and  donor  stack.  Figure 
9  below  shows  how  this  is  done. 


Acceptor 


Ds  =  diameter  of  shield  (bar)  or  width  of  shield  (plate) 
t  =  thickness  of  shield 

=  distance  of  shield  from  donor  source  (as  shown) 

$2  =  distance  of  acceptor  group  from  shield  (as  shown) 

($1  +  t  +  S2)/(Si  +  t/2)  =  TWDTH/Ds 

Si(TWDTH  +  Ds)=  Ds(S2  +  t)  -  t/2(TWDTH)  (25) 
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Step  4» 


On  a  scaled  drawing  (as  shown  in  figure  A.l)  with  donor  and 
acceptor  stacks  spaced  the  distance  S  and  their  bases  on  the 
same  horizontal  plane,  lay  off  the  dividing  rays  for  polar 
zones  of  10°  centered  at  integral  multiples  of  10°. 
Determine  zones  in  which  fragments  can  hit  the  acceptor 
projectiles,  ignoring  the  base  plate  and  fuze  section  if  no 
funnels  are  involved.  As  stated  earlier,  estimate  the 
fraction  F2i  of  the  zone  over  which  the  vulnerable  part  of 
the  acceptor  extends.  F^i  =  1.0  for  zones  completely 

spanned.  If  S  is  large  enough  that  only  the  central  zone  is 
involved,  then 


Fz  =  ISH/irS 


where  H  =  height  of  the  vulnerable  portion  of  the 
projectile. 

If  funnels  are  involved,  then  H  is  the  combined  height  of 
the  projectile  and  funnel,  ignoring  only  the  base  plate. 

Also  estimate  the  average  outside  diameter  of  projectile 
(in),  casing  thickness  (in),  and  the  angle  of  impact  at  the 
surface  of  the  shield  for  each  zone. 


Step  5. 


Assuming  that  there  is  no  decrease  in  the  weight  of  a 
fragment  after  perforation,  the  iterative  process  is  greatly 
simplified.  For  each  zone,  compute  the  critical  mass  mf.-j 
(grain)  by  iteration  from  the  following  equations: 


=  e(-0.0025S2mi-l/3) 

=  e(-0.0025Simi-l/3) 

ViNi  =  VoFv-j 
Vspi  = 

A  =  (m-j/k)2/3 

Vp.  =  10^1(tA)“lmi^l(sec  6)^1 
Vri  =  Vspi  -  10C(tA)“mi^(sec  e)^V^p. 
VsH  "  VriFoi 


(20) 

(26) 

(27) 

(28) 


(29) 

(30) 

(31) 

(32) 
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(22) 

(13) 


d  =  (4F/\/Tr)l/2(ni-j/l<)l/3 
VsH  =  (1  +  Kix/d)Jd-a 

Starting  with  m-j  =  itIq,  if  Vp.  >  Vgp.  for  any  value  of  m-j , 
then  the  fragment  does  not  perforate  tlie  shield.  By  setting 
m-j  to  be  equal  to  a  large  value  (5,000  grains),  the 
probability  of  such  a  fragment  with  zero  residual  velocity 
causing  detonation  is  automatically  set  to  zero. 

Step  6. 

For  each  zone  also,  compute  the  total  n^^mber  of  fragments 
N-j ,  the  azimuthal  factor  F(^-j  and  fraction  f-j  of  the  striking 
fragments  that  have  mass  exceeding  the  critical  mass  m-j 
using  these  equations: 

Ni  =  CAFm(0)mo  (17) 

Fai  =  2(D  -  x)/2  (Sj  +  e  +  S2)  (33) 

f.  =  Fie-(2niri/nio)^/^  (16) 


Step  7. 

The  probability  of  detonation  is  calculated  as  follows: 


N  =  ^  NiFctiFzifi  (14) 

i=l 

P  =  1  -  e-N  (19) 

The  notations  above  have  been  defined  previously  in  an 
earlier  section  of  the  report. 
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DISCUSSION  OF  RESULTS 


Although  the  second  part  of  the  report  dealt  mainly  with 
stacked  projectiles,  calculations  were  done  to  predict  the 
probability  of  detonation  propagating  between  adjoining 
individual  155mm  (M107  Comp  B)  projectiles  with  and  without  plate 
shields  between  the  projectiles. 

With  no  shield  present,  the  model  predicted  a  probability  of 
detonation  of  49  percent  for  a  distance  of  0.61  m  (24.1  in),  but 
the  confidence  level  could  not  be  determined.  However,  when  a 
0.5-inch  thick  steel  plate  was  placed  between  the  two 
projectiles,  the  probability  of  detonation  was  zero  for  the  same 
distance.  Tests  were  done  with  the  same  projectiles,  same 
thickness  of  shield  and  same  spacing  between  projectiles,  and  out 
of  48  trials,  no  detonation  of  the  acceptor  projectile  was 
observed  (ref.  14).  The  number  of  tests  is  enough  to  give  a  high 
confidence  level.  Sample  calculations  are  presented  in  appendix 
A. 


Calculations  were  also  done  for  groups  of  81mm  and  105mrn 
projectiles  to  determine  the  adequacy  of  the  model.  No  favorable 
correlations  were  expected  since  parameters  designed  for  single 
and  groups  of  four  (2x2)  projectiles  were  used  in  the  model  to 
predict  the  probability  of  explosion  propagation  between  groups 
of  16  (4  X  4)  lU5mm  and  72  (6  x  12)  81mm  projectiles. 

For  a  separation  distance  of  109. 7  m  (360.0  in)  between  two 
groups  of  16  (4  x  4)  105mm  (Ml)  projectiles,  the  model  predicted 
a  probability  of  detonation  of  20.6  percent  compared  to  10 
percent  at  a  95  percent  confidence  level  provided  by  test 
results.  Similarly,  the  model  predicted  a  probability  of 
detonation  of  2.4  percent  for  a  separation  distance  of  0.914  m 
(36.0  inches)  between  two  groups  of  72  (6  x  12)  81mm  (M374A1) 
projectiles.  Test  results  gave  a  probability  of  6.8  percent  at  a 
95  percent  confidence  level  (refs.  13  and  15). 

Finally,  the  model  was  used  to  predict  the  probability  of 
detonation  for  individual  8-in  (M106)  HE  projectiles  spaced  at  a 
distance  of  0.305  m  (12  inches)  with  a  3-inch  diameter  aluminum 
bar  placed  between  the  donor  and  acceptor  projectiles.  The  model 
came  up  with  a  probability  of  zero.  Test  results  show  that  out 
of  50  observations,  no  propagation  was  observed  thus  providing  an 
upper  limit  of  7.1  percent  for  a  confidence  level  of  95  percent 
(ref.  16). 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Concl usi ons 

The  model  and  procedures  described  in  this  report  provide  a 
rational  basis  for  determining  the  probability  of  detonation 
propagation  between  adjoining  explosive  items.  From  the  analysis 
of  available  data  and  calculations  provided  in  this  report,  the 
following  conclusions  have  been  reached: 

1.  When  used  in  a  generalized  way,  the  model  can  be  a 
valuable  tool  in  the  design  and  layout  of  facilities 
for  the  manufacture  and  storage  of  explosive  items. 

2.  To  compute  the  required  separation  distance  S  between 
projectiles  for  a  given  probability,  the  model  can  be 
used  to  determine  the  probability  of  detonation  for 
trial  values  of  S  until  a  correct  value  is  found.  This 
was  determined  to  be  the  most  accurate  means. 

3.  Certain  approximations  and  assumptions  made  affect  the 

accuracy  of  the  model,  especially  when  stacked 

projectiles  are  involved.  The  accuracy  would  be 
enhanced  by  appropriate  experiments  such  as  those 
recommended  below. 

Recommendations 

The  following  recommendations  are  made: 

1.  Perform  further  analyses  of  this  modified  model  in  the 

areas  that  cover  shields,  funnels  and  stacked 

project! 1 es. 

2.  Design  a  test  plan  covering  areas  which  should  be 
further  investigated  experimentally  such  as  arena, 
fragmentation  and  sensitivity  tests. 

3.  Perform  arena  tests  for  projectiles  stacked  in 
configurations  of  interest,  with  the  view  of  collecting 
information  on  velocity,  polar  and  azimuthal 
distribution  of  the  fragments. 

4.  Perform  fragmentation  tests  to  determine  the  number  of 
fragments  emitted  from  different  configurations  of 
stacked  projectiles  (and  the  distribution  of  these 
fragments) . 
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5.  Perform  sensitivity  tests  (similar  to  the  ones  begun  in 
ref.  12)  to  accurately  determine  the  relationship 
between  boundary  velocity  and  each  of  the  following: 

a.  Fragment  mass 

b.  Presented  (or  impact)  area  of  fragment 

c.  Thickness  of  acceptor  plate. 
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APPENDIX  A 
EXAMPLE  PROBLEMS 


Sample  Calculations 

1.  81mm  (M374A1)  Projectile  with  No  Shields 


Average  Outside  Diameter 
Average  Casing  Thickness 
Height  of  Charge  in  Projectile 
Density  of  Enclosed  Charge  (Comp  B) 
Density  of  Metal  (Steel) 


D-j  =  3.04  inches 
X-j  =  0.29  inch 
H  =  7.2  inches 


Pc  =  0.058  lb/in3 
Pm  =  0.283  lb/in3 


Mass 

of  metal 

(Sidewalls)  M 

=  5.0  lb 

=  35,000.0 

grai ns 

Input 

Table  (A.l). 

81mm  Projectiles 

Pol  ar 

Casi ng 

Outside 

angl  e 

thickness 

di ameter 

Zone 

 Jdej) 

FZ 

FV 

FM 

( 

1 

80 

0.57 

0.862 

0.183 

0.29 

1.91 

2 

90 

1.00 

0.967 

0.377 

0.29 

3.04 

3 

100 

0.43 

0.969 

0.377 

0.29 

3.08 

Note:  The  casing  thickness  and  outside  diameter  for  each  zone 
are  measured  at  the  point  where  the  centerline  of  that 
particular  polar  zone  intersects  the  centerline  of  the 
projectile. 
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FIGURE  A.I  81  mm  (M374AI)  projectiles 


Input  Format-  Example  1 
Card  Type  1 


NCASE  =  0 
NPROJ  =  1 

Card  Type  2 

Title:  81-mm  (M374A1)  Projectiles  (Comp  B) 

Card  Type  3 

Type  of  Shield:  None 

Card  Type  4  (Needs  decimal  points) 

VC  =  8800.0 
DENC  =  0.058 
DENM  =  0.283 


SK 

=  660.0 

XM 

=  35000.0 

H 

=  7.2 

XD 

=  0.29 

DD 

=  3.0 

Card  Type  5 

CO 

=  4500.0 

SB 

=  2140.0 

Cl 

=  0.8 

FA 

=  0.5 
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FI  =  0.2 
PLATE  =  0. 

Card  Type  6  (A  blank  card  can  be  substituted) 

ALPHA  =  0 
BETTA  =  0 
GAMMA  =  0 
CONST  =  0 
XLAMD  =  0 

Card  Type  7  (A  blank  card  can  be  substituted) 

ALPHl  =  0 
BETAl  =  0 
GAMAl  =  0 
CONTI  =  0 

Card  Type  8 

SA  =  0.61 
CB  =  0.283 
XJ  =  1820.74 

Card  Type  9 
S  =  21.08 
NZONE  =  3 
NMIS  =  1 

TWDTH  =  blank.  This  value  can  be  left  out  since  it  only 
applies  to  stacked  projectiles. 
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Card  Type  10  (Values  are  taken  from  Input  Table  A.l).  One 
card  required  for  each  zone. 


1st  Zone 

2nd  Zone 

3rd  Zone 

FZ  =  0.57 

FZ  =  1.0 

FZ  =  0.43 

FV  =  0.862 

FV  =  0.967 

FV  =  0.967 

FM  =  0.183 

FM  =  0.377 

FM  =  0.377 

XR  =  0.29 

XR  =  0.29 

XR  =  0.29 

DR  =  1.91 

DR  =  3.04 

DR  =  3.08 

ETTA  =  1.0 

ETTA  =  1.0 

ETTA  =  1.0 

THETA  =  - 

THETA  =  - 

THETA  =  - 

Card  Type  11  (End  of  Data  Indicator) 

A  blank  card. 

2.  155mm  (M107)  Projectiles  with  0.5-Inch  Steel  Plate 

Shield. 


Average  Outside  Diameter 
Average  Casing  Thickness 
Height  of  Charge  in  Projectile 
Density  of  Change  in  Projectile  (Comp  B) 


D-j  =  6.0  inches 
X-j  =  0.29  inch 
H  =  7.20  inches 
P  =  0.058  lb/in3 
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Density  of  Metal  (Steel) 
Mass  of  Metal  (Sidewalls) 


p  =  0.283  lb/in3 
77.5  lb  =  542,500.0  grains 


Input  Table  (A. 2).  155mm  Projectiles  with  shielding 

Pol  ar 


Zone 

angle 

(deg) 

FZ 

FV 

FM 

XR 

DR 

n 

6 

(deg) 

6 

(rad) 

1 

70 

0.07 

0.732 

0.07 

0.68 

4.09 

1.0 

20 

0.35 

2 

80 

1.00 

0.862 

0.18 

0.68 

5.45 

1.0 

10 

0.17 

3 

90 

1.00 

0.967 

0.377 

0.68 

6.09 

1.0 

0 

0 

4 

100 

1.00 

0.967 

0.377 

0.82 

6.00 

1.0 

10 

0.17 

5 

no 

0.87 

0.862 

0.18 

1.36 

6.54 

1.0 

20 

0.35 

Input  Format  -  Example  2 
Card  Type  1 

NCASE  =  0 
NPROJ  =  1 

Card  Type  2 

Title:  155mm  (M107)  Projectiles  with  0.5-inch  Steel  Plate 
Shield 

Card  Type  3 

Type  of  Shield:  0.5-inch  Steel  Plate 

Card  Type  4 
VC  =  8800.00 
DENC  =  0.058 
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24,1 


.0>2 
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FIGURE  A.2  I554nm  (MI07)  projtctites  with  0.5-  inch  steel  plote. 


DENM  =  0.283 


SK  =  660.0 
XM  =  542500.13 
XD  =  0.68 
DO  =  6.09 

Card  Type  5 

CO  =  4500 
SB  =  2140 
Cl  =  0.8 
FA  =  0.5 
FI  =  0.2 
PLATE  =0.5 

Card  Type  6  (From  Table  1) 

ALPHA  =  0.889 
BETTA  =  0.945 
GAMMA  =  1.262 
CONST  =  6.390 
LAMDA  =  0.019 

(Beta  values  are  entered  as  positive  values  always) 

Card  Type  7  (From  Table  2) 

ALPHl  =  0.906 
BETAl  =  0.963 
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GAMAl  =  1.280 


CONTI  =  6.52 

Card  Type  8 

SA  =  0.61 
CB  =  0.283 
XJ  =  1820.74 

Card  Type  9 

51  =  7.5 

52  =  16.1 
NZONE  =  5 

NMIS  =  1 
TWDTH  =  (blank) 


Card  Type  10  (From 
zone* 

Input  Table  A.2) . 

One  card  for  each 

1st  Zone 

2nd  Zone 

3rd  Zone 

FZ  =  0.07 

FZ  =  1.0 

FZ  =  1.0 

FV  =  0.732 

FV  =  0.862 

FZ  =  0.967 

FM  =  0.07 

FM  =  0.18 

FM  =  0.377 

XR  =  0.68 

XR  =  0.68 

XR  =  0.68 

OR  =  4.09 

DR  =  5.45 

DR  =  6.09 

ETTA  =  1.0 

ETTA  =  1.0 

ETTA  =  1.0 

THETA  =  0.35 

THETA  =  0.17 

THETA  =  0 

4th  Zone 

5th  Zone 

FZ  =  1.00 

FZ  =  0.87 

FV  =  0.967 

FV  =  0.862 

FM  =  0.377 

FM  =  0.18 

XR  =  0.82 

XR  =  1.36 

DR  =  6.0 

DR  =  6.54 

ETTA  =  1.0 

ETTA  =  1.0 

THETA  =  0.17 

THETA  =0.35 

Card  Type  11  (End  of  Data  Indicator) 
A  blank  card. 
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3. 


lOSmm  (Ml)  Projectiles  -  Sixteen  (4  x  4)  on  Pallets 


I®. 


Average  Outside  Diameter 
Average  Casing  Thickness 
Height  of  Charge  in  Projectile 
Density  of  Charge  (Comp  B) 
Density  of  Metal  (Steel) 

Mass  of  Metal  (Sidewalls) 


D-j  =  4.133  inches 
Xi  =  0.489  inch 
H  =  11.8  inches 
Pc  =  0.058  lb/in3 
pfji  =  0.283  Ib/in^ 

M  =  25.5  lb  =  178,500.0  grains 


Determination  of  C/\ 
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For  an  azimuthal  zone  bounded  by  0  and  180°  line,  the  shaded 
regions  show  the  portions  of  projectiles  that  are  not  involved  in 
the  fragment  mass  distribution.  The  equivalent  weight  of  5 
projectiles  is  involved  in  fragment  mass  distribution: 

C/\  =  5  X  Wt.  of  single  105mm  projectile  (grains) 

=  5  X  178,500.0  grains 
=  892,500.0  grains. 

Determination  of  FZ: 

For  a  separation  of  360  inches,  only  the  middle  (90-degree) 
zone  is  involved. 

FZ  =  18H/TTS  =  (18  X  11.8)/tt  x  360  =  0.188 


Input  Table 


Table  A. 3.  105mm  (Ml)  Projectiles 


Pol  ar 
angl  e 


XR  DR 
(in)  (in)  n 


Zone  (deg)  FZ_ R 


FM 


1 


90  0.188  0.886  0.206 


0.489  4.133  2.0 


Separation  distance  =  360  inches 
(See  notes  at  bottom  of  Table  A.l) 
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Input  Format  -  Example  3 
Card  Type  1 


NCASE  =  0 
NPROJ  =  1 

Card  Type  2 

Title:  105mm  (Ml)  Projectiles  Grouped  16  (4  x  4)  on  Pallet 

Card  Type  3 

Type  of  Shield:  None 

Card  Type  4 
VC  =  8800.00 
DENC  =  0.058 
DENM  =  0.283 
SK  =  660.0 
XM  =  892500.0 
H  =  11.8 
XD  =  0.489 
DO  =  4.133 

Card  Type  5 
CO  =  4500.0 
SB  =  2140.0 
Cl  =  0.8 
FA  =  0.5 
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FI  =  0.2 


PLATE  =  0.0 

Card  Type  6 

ALPHA  =  0 
BETTA  =  0 
GAMMA  =  0 
CONST  =  0 
LAMDA  =  0 

Card  Type  7 
ALPHl  =  0.0 
BETAl  =  0.0 
GAMAl  =0.0 
CONTI  =  0.0 

Card  Type  8 
SA  =  0.61 
CB  =  0.283 
XJ  =  1820.74 

Card  Type  9 
S  =  360.0 
NZONE  =  1 
NMIS  =  4 
TWDTH  =  19.53 
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Card  Type  10  (From  Input  Table  A. 3).  One  card  required  for 
each  zone. 

1st  Zone 
FZ  =  0.188 
FV  =  0.886 
FM  =  0.206 
XR  =  0.489 
DR  =  4.133 
ETTA  =  2.0 

Card  Type  11  (End  of  Data  Indicator) 

A  blank  card. 
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APPENDIX  B 

INPUT  FORMAT  FOR  COMPUTER  PROGRAM 


Introduction 

Tfiis  section  presents  the  formats  used  for  specifying  the 
various  input  parameters.  Each  type  of  card  is  described  below 
in  terms  of  data  format,  definition  and  field  locations.  The 
numbers  above  the  graphic  representation  of  each  card  identify 
the  last  column  in  each  data  field  of  that  card.  The  letters 
below  designate  the  format  for  the  data.  In  fields  designated 
"I",  the  quantity  entered  must  be  right-adjusted  to  the  last 
column  in  the  field  and  cannot  contain  any  decimal  point.  In  the 
fields  designated  "F",  a  decimal  point  is  required;  however,  the 
number  can  be  located  anywhere  in  the  field.  Fields  designated 
"A"  are  alphanumeric  fields  and  have  to  be  left-adjusted, 
beginning  at  the  first  column  of  the  field. 

Data  Cards 

Card  Type  1 


=  1  if  perforation  of  the  projectile  casing  is  the  failure 
criterion 

NPROJ  =  number  of  times  the  program  is  to  be  used.  If,  for 
example,  calculations  are  required  for  the  01mm  and 
155mm  projectiles,  then  NPROJ  =  2.  Note  that  NPROJ  does 
not  mean  the  number  of  projectiles  in  a  stack. 

Card  Type  2 


TITLE 


"A"  Format 
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Card  Type  3 


16 

SHIELD 

■A'  FORMAT 

If  no  shields  are  present,  write  "None".  If  this  card  is  left 
out,  when  no  shields  are  present,  erroneous  results  will  be 
obtai ned. 


VC 

DENC 

DENM 

SK 

XM 

H 

XD  :  DD 

_  M 

F"  Format - » 

VC  =  Gurney's  constant  (fps) 

DENC  =  Density  of  charge  (Ib/in^) 

DENM  =  Density  of  metal  (Ib/in^) 

SK  =  Shape  factor  (grain/in^) 

XM  =  Weight  of  metal  (grain) 

H  =  Height  of  charge  in  projectile  (in) 

XD  =  Average  casing  thickness  (in)  measured  at 
middle  of  projectile 

DD  =  Outside  diameter  of  projectile  (in)  measured 
at  middle  of  projectile  also. 

Card  Type  5  ' 


10,  20.  30,  40,  50  60 


CO 

SB 

Cl 

FA 

FI 

PLATE 

^ - iirii  _ .4-  - 

"F"  Format 


CO  =  Average  mass  coefficient 
SB  =  Sensitivity  parameter 
Cl  =  Protection  coefficient 
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FA  =  Area  factor 


FI  =  Impact  factor 
PLATE  =  Thickness  of  shield  (in) 


Card  Type  6 


10 

20 

30 

40 

50. 

ALPHA 

BETA 

GAMA 

CONST 

XLAMD 

y/////////A 

"F' 

"  Format 

These 

are  constants  obtained  from  table  1  of  the  report. 

Card  Typ 

e  7 

20 

30 

40 

ALPHA  1 

BETA  1 

GAMA  1 

CONT  1 

y/////////////. 

Constants  taken  from 

^ - 

table  2. 

"f"  Format 

When  no  shields  are  present,  the  constants  of  card  types  6  and  7 
should  contain  zeros.  An  erroneous  result  will  be  obtained  if 
these  cards  are  left  out  entirely. 


Card  Type  8 

10  20 


SA 

CB 

X  J 

30 


Y//////////////Z\ 


Constants  used  in  equations  9  and  13. 
CB  =  B  in  equation  9. 


T"  Format 


SA  =  a) 

)  equation  13 
XJ  =  J) 


If  no  shields  are  present,  then  card  type  9  would  be 


10  20  30  40 


S 

NZONE 

NMIS 

TWDTH 

Y////////////, 

h - -'I" - -pF"— ^ 

If  shields  exist  between  donor  and  acceptor  projectiles,  then 
card  type  9  becomes 


SI 

S2 

NZONE 

NMIS 

TWDTH 

y///777//. 

- ^^'F" - -j 

- - "I' - 
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S  =  Center-to-center  spacing  between  donor  and  acceptor 
projectiles  (in)  for  no  shields  present. 

51  =  distance  from  centerline  of  donor  projectile  to  inner  face 

of  shield  (in) 

52  =  distance  from  outer  face  of  plate  to  centerline  of  acceptor 

projectiles. 

NZONE  =  Number  of  zones  involved  (see  Sample  Calculations). 

NMIS  =  Number  of  projectiles  in  a  row  or  column  facing  the 
donor  source. 

TWDTH  =  Length  of  the  row  or  column  of  projectiles  facing  the 
donor  source. 


Card  Type  10  (one  required  for  each  zone) 


10 

20 

30 

40 

50 

60 

70 

FZ 

FV 

FM 

XR 

DR 

ETA 

THETA 

"F"  Format 


FZ,  FV  and  FM  are  parameters  defined  in  the  report.  XR  and 
DR  are  the  casing  thickness  and  outside  diameter  in  inches 
of  acceptor  projectile.  The  value  of  ETTA  can  be  obtained 
from  the  second  part  of  the  report.  THETA  is  the  fragment 
impact  angle  of  the  face  of  the  shield. 

End  of  Data  Indicator 


A  blank  card  is  required  at  the  end  of  data. 

Program  Listing 

The  PODS  Program  is  written  in  FORTRAN  language  for  the  IBM 
1130  Computer.  Modifying  it  to  run  on  any  system  should  not  be  a 
difficult  task. 

The  program  is  structured  to  predict  the  probability  of 
detonation  using  either  of  the  following: 

1.  Perforation  of  projectile  casing  as  a  failure  criterion 
(NCASE  >  0). 

2.  High  order  detonation  as  failure  criterion  (NCASE  =  0). 
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Because  of  the  dimensional  constraints,  only  twenty 
10-degree  zones  can  be  involved.  This  imposes  a  restriction  to 
how  close  one  projectile  can  be  to  another;  however,  the 
distances  involved  in  ammunition  facilities  are  within  the 
1 imits. 

As  explained  under  "DESCRIPTION  OF  THE  MODEL"  in  the  main  body 
of  the  report  (on  page  4),  the  computer  program  is  presented  below. 
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IF  (NCASE-I ) 103, 104,1 OA 
103  CALL  A0LRO(FINAL,SX) 

DI ST=SX 


o  ♦ 

1^  • 


o  ♦ 

• 


o  ♦ 

fO  • 


o  ★ 

rj  • 


00 


o 

oc 

0- 


c 

f\J 


o 

c\i 


< 

I— 

00 


Z  Qi  f\4 
l-H  O  t- 


KX  ^ 

I—  N-  00 
00  t-  *-4 


c 

fM 


^  O 
X  O 


M-(  ^  «r  t- 


1-4  1-4  ^ 

O  O  X  z 

o  o 

O  O  <NJ  >-H 

•  o  o  ^>— 

•  o  c  -  <r 

♦  fO  O  Z  » 

•  X  V  K  -  O  LU  » 

•  00  w  »—  Ui  Z  I 

•  II  LuOLU^C^i-H 

^  •00t-(  I-HQCU-Z. 

•  I-HQCOQCOOO 

•  Q  r*  o  *:»  u.  tj  « 

■n  ^  o  ^ 

•  4-0 


CO  O 
O  <M 
a  W 

a.  > 

X  u. 

w  K 
Q. 

a  O 
-J  (\i 
O 

<  rJ 


z  z 

l-H  O 
»—  >-H 

O  00 

o  z 

0£  UJ 


<  Ui  <  < 
u.  ^  o  O 

V  <  \  s 

X  -J  <r  T- 
00  a.  I—  <r 

V  ^  I—  I-  <c 
ffi  O  Ui  LU  £ 

00  z 

^  X 

O  V 

u>  C3J 


z 

X 


4-  O 

<-i  > 

z  z 
o  o 
z  t: 
^  s: 
o  o 
u»  u» 


K  V  V 

<  r- 
XXX 

a.  a.  ^ 

«|  _|  <  vO 
<t  <  00  4- 

Z  Z  Z  4-  ' 
O  O  O  •  I 
»:  ar  z  fo  I 

z  ar  z:  n  i 

0  0  0  1-41 

o  (->(->  a.  : 


X 
/-s  in 
O  4- 

O  \ 

o  - 
>t  o 

\  ■ 


^  -tt 
-  41 

^  4t 
00  « 

I  z  ♦ 

.  I-H  4t 
4t 

QC  ★ 
UJ  -H 

I—  ★ 

UJ  4t 

ac  ♦ 

j  <  ★ 
:  i-H  -tt 

>  o  ★ 
J  -n 
QC  -li 
UJ  -d 
»-  -d 
O  -d 
O  4t 


in  ' 


00 

z 

o 


O  UI  z 
Z  -J  o 

O  i-H  ►H 
u  I—  I— 
u>  < 
00  UJ  z 
<t  ->  o 
O  t- 
Z  QC  UI 

o  a.  o 

1-4 

»—  UJ  CX 
<  X  Ui 


»—  u.  O 
UI  o 
O  X 

z  o 
u.  O  »-i 

O  i-i  X 
I—  • 

>-  <t  z 
»—  ex  o 

1-4  O  »-4 
^  u.  q: 

»-4  oe  UI 
CD  UJ  >— 

«t  Q_  1-4 
X  OC 

o  ^  u> 
ex  tx  o 

OL  UJ  Ui  Ui 

>  q:  Qc 

UJ  Ui  o  UI 

X  **  —I  o 


-d 

00  -d 

:  z  41 

I-H  4t 
1  W  4i 
00  4c 


.  O  l-H 
X  <  * 


00 


z 

o 


I—  •  < 

<  o 

-J  z  o  z 

3  UJ  O 
UJ  O  ex  M 
^  ^  Ui  qe: 

<  Ui  o  Ui 
UJ  i-H  »-4  ^— 

X  00  1-4 

00  z  ex 
a,  o  cj 


o  <  1“  e 
<  o  ; 
Ui  z  . 

Ui  UJ  I 

z  z  00  . 

1-4  <  M  I 
I—  I— 

3  t/J  O  I 
O  *-4  Z 
ex  ea  i-H  I 
eu  t/3 

3  u.  «t 
41  00  O  UJ 


i-H  in  ^ 


Z  O 

K4-  > 


«  (\J  Ui  u. 


WZ 
4-  O 
\  rsj 
4-  Z 
N 

O  00 


1-4  4(  4t  4t  i-H  ex 


1-4  Z  Z 

CJ  O  »si  ^  o 

\  fsi  u.  in  rsj 

in  z  ^  •  z 

•  K4-  o  ^ 

0  4-  0  4-  4- 

4-  II  o  u.  o  II 

U-  l-H  fO  00  •  _J 

,u^  4.W  O  O 

»-  N.  rj  I-  II  4-z- 

<  ^  w  <  ex  c  X 

at  r-  o  ar  <c  oj  ii 
ex  <  ex  CD  4- 

O  O  UJ  O  Z  O  S' 

u.  c:^  ex  u.  X  o  X 


KO  > 

ro 


•n 

KO 


o  ^  ' 


X  4t 


«l  41  4t 


4i  I—  4t  «r 


\  00 
-J  4t 
Z  O  l/^ 
^  O  <M 

4-  OJ  O 

O  K  o 

0  4-  • 

^  II  c 
4.x  I 
in 

m  Q. 
UJ  c  X 
^  I—  Od  UJ 
tl  1-4  II 
— I  CX  O 


I-  1^ 
Ui  <M 
4t 


1^ 

O 


kO 

KO  4- 

O  4- 

r\j  K 

oc 

4-  O 


4t  <\J  I 
^  41  I 
-J  -d  ' 
X  ' 
>  to 

u.  W 

4t  w 
O  I 
7>  4-  I 
II  II 


7  3  Xi  U-  >  >-  X 


O 

O 

O 


4-  O 

O 

O 


N-  4- 

^  O 

1^  o 

KO 


y  1-4 

X  *—  o 

i  W  o 

|4"|  oO  4-  K' 

X  I  z 
I  X  <r  b-e  X 

II  w  w  II 
I  X  4- 

t-4  u-  u.  y 
^  »-l  X 

r-  X 
c  cr 


UJUJUJUJUJUJUJUJUJUJ 


o 

UJ 


61 


o 

o 


2 

Z) 

o 


r-OJrO'>-Tir\'Or^cx30o^f\iro>4’U^>Oh^oooat“r\jfO'«4’u^'Of^of>oat-r\jrO'«4’u^'Or^oooot“fMrO’^u>'Or^oorK 

Lnu>ir\mir\ir\ir\u>u>>0'0'0'0'0'0'0'0'0'0r*-r^r^f^h-f^f^f^f^r^ooooooooooooooooooooo'0'0'0'0'0os0'0o» 


o 

o 


^  ♦  Li_ 


t-t  •'N.  tn  ^ 


K',  ^ 

ro 

V  _l 
O  ^ 
ro  Qi 
ro  X 
o  \  -tt 

r*  O  rsi 
X  m  I 

'v.  f«r«  ^ 

o  ^  ^  _i 

a  _l  r-  W 

«“  r-  ^  I  a  (\i 

^tr  ^OFTi/'oro 

5  rvj  i/^  rvi  u.  1-1  ' 


s:  I 

I  X  , 


2 

I  I 


CD 

2 


-  II  .. 

i-joujor:2  •to<^-»Qc 


.  h-  u-  U-  <  I— 


I  o  o:  o  2  -T  o 

I  O  *!1  X  I~I  U.  O 


u.  (/)  X  O 


2  ^ 

\Ln  \ 

LU  •  o. 

2  0  t-H 

O  ^ 

NJ  U.  > 

SC  V  U. 

SCD  \ 

f\j  «“ 

to  *-H  Ui  *-H 

\  f\J  2 

r-  V  LU  O  M  <-N 

I  CO  in  h-  LI.  i/> 

^  •  -r  ar  ^  • 

Ol  fvj  O  _i  Kro  O 

X  O  «“  Q_  r-  O  «“ 

lUf\JOLi.-f  IIOU-O 

irororvjfvjp-irooo  • 

r-fVJ  K'-'CO  \^0 
y5ll  fVJI—  II 

II  a:o^-rr-c:^<ca 
Oh- oart/)r^os:< 
a  <*  a:  II  a:  CD 

CL  OuJOh-OLUOS: 
xOofu.  cooca^Li.x 


ro  fO 
ro  ro 
ro  ro 
ro  fsn 

0^0 


f^ 

nO 

o 

o 


r-  OX 


»< 
I  2 


f\J  «“ 

to  to 
*  * 
to  to 
a  ro  f\j 

a  a  a 

lo  CD  a 


a  >—  c  <  <  o  r- 

♦  2  Q.  I—  I-  -It  s 

■K  O  Cl  LU  LU  X  >0 
^0<CDX2_I«- 
^♦«_»-fch-Z)LU»-» 
to-»t  -It  ♦'“-'cn^a 


a 

a 


a  a 

I  I 


I  O  r- 

tj  z: 


^  Ui  r-  to  -It  Ol 
O  h-  Z"  O  X  X  ( 


II 


W  w  -n  >- 


lO  o 


O  _J 
O  2 


Q-  d.  I 

X  X  I 

LU  X  r^  LU  Ui  I 

I-  II  r-  V  II  11 
*-H  «“  II  O  O 
or  y  o  ^  o  o  I 

3  >-  cj  ll.  o  : 


2 

r-t  II 
I  >  < 
II  Q- 
Q_  Q_ 
I  to  < 

>  O 


r-  r-  r-  2C  to  u. 
^  ^  w  II  >  ►_! 

II  II 
X  X 


t-i  13 

CT  CD 


II  II  _l  II  ^ 

X  X  LU  U.  w 

2  LD  >  U. 

-r  2  Ol  *-H  u- 

o  <  X  £3  r-i  >- 


s: 

Z" 

< 


^  ♦ 
to  <r  h-  <  2 

2  CL  h-  h-  O  < 

O  OL  LU  UJ  Z~  O 

O  <t  2D  X  <t  ♦ 

■H  tJ  -k  I—  _J  2  LU 

•  ^t-k-k'^xibz: 

a  <-»Uiit-tO-k2D<0 

O'©  oi-z’o-k-kh- u- 
OCIU  •-r>-oCL2l  -k 
CDlMDCD—l'^'^tOr-iOL'-H 
O  2»-CLr-r->02loa: 
t-Or-*'-''-''-''-'!!  II  >> 
II  I—  h-  II  II  II  II  O  LU  II  II 

r-  2222tO_iart-<X 
Z"OOr-«I3<2LU<'tytO 
t5oancr^<>>—  >> 


ro  o  ro 

o  r-  c 

ro  r-  ro 


C  r-  ro  «- 

ro  r<^  ro  o 

ro  ro  ro  ro 


ro  ro 
a 
o 


ro 
o  a 
1^  o 
ro 


oc 


2 

Z" 

D 


O 

LJ 


62 


2 

z: 

Z3 


O 

O 


3 

O 

o 


T-rjro>j'm>or^cooat-ourO'^iA'«or«»oooaT-rjro>j-m>or^coooT-rjro>j-iA'OrwaoooT-f\jf>o>a-iA>oN.jooo 

oooaaoaoo<-t-T-T-T-T-t-t-T-T-f\jrarjrjrjrarjrjrjfNjrororofOrorororororo>3->3->3->j->3->3--st>t>t>3-iA 


o  ♦ 

00 


c  ♦ 

ru 


ro 
fo 
i/>  ro 
•  ro 

o  • 

♦  o 

♦  ♦  ^  ♦ 

1-4 

^  O  LU 
Q_  w  >— 


lO 

sr  o 
c  ♦  I 

HH  : 

o  *:  : 


♦  y* 
>5'  >- 


I  o  cr  I 

.  ♦  X  T-  : 

•  y  ♦  w 

c  t-  ♦  ' 

LU  O  00  ' 

II  II  00 
y  LU  II  I 

^  ^  ^  I 


«t  «t  I 
UJ  ^  I 


rj 

•  o 

O  t- 
lO  MO 

I  nT  ■ 

cr  <->  V  t— 

oo  I—  ^  t- 

>  o  ^ 

I  — I  O  X 

I—  00  o  -> 

q:  o  +  I 

00  00  t-  o  : 

00  I  y  : 

II  «x  X  X  I 

I-  N-'  W  tl  I 

o  T-  : 

_i  u-  u_  y  o 

on  HH  HH  >-  o 

ro  >«•  N- 


0- 

♦ 

ru 


ro  /-> 
ro 

V  _J 

to 

ro  oc 
ro  X 
o  V  ♦ 
y  ro  fM 
X  ro  I 
ro  ^ 

o  ^ 

O  —It-'*-' 

*o  T-  rvi  N-  I  Of  lo 

»  o  V  fM  y  00  o  ro 

>  fM  in  fM  U_  K-i  w  K-( 

:  ♦  y  II 

I  o  UJ  o  y  2  X  o 

'  I—  I—  X  w  u.  »— 

:  hh  II  _i 

>  o  a  o  >-  u-  c  o 

»  09  3  0^  X  HH  U.  UJ 


w  «  I 
Q1  /-N 

X  o  : 


y  >-  ; 
z  *  ; 
I  fM  > 
3:  '*-' 

r-  I  i 

O  *«-o  < 

3  Ol  I 


.  C 

Ll. 

'  -n 


ro 

o 


X 

UJ  r*- 

09  T- 

q:  \ 

UJ  ro 

>  • 

Z  CD 
O  t- 
O  u- 


LU  >0 
O  U_ 

O  V 


UJ  ro 

o  ■ 

O  >0 
Q1  U. 
0-  V 
X 
z  >«• 

O  V 

HH  ro 

t—  • 

.  <  NO 


II  « 
X  I-I  ( 


00  I-  X 

y  HH  >«• 

X  2  V 
1  lO  fM 

**-'  ro  »-H 

a.  V  V 

X  -  X 

UJ  t- 

I  -  t-  I 

I  X  I—  T-  '*-'  N-'  : 


C2 

(M 


\  t-i 
u. 
O  \ 


O  T- 

y  I- 

c  z 

-I  o 

X  o 

s  \ 

y  c  T- 
X  y 
V  y  y 


CE 

o 

O'  *>^  U.  r—  09  09 
Q_  >  \  C  \  \ 

rsj  u_  ^  ^  C  T- 

\  00  Q_  r—  < 


-I  II 
C  00 


o 

o 


NO 

o 

fM 


u-  u-  X  O 
>«•  to  T- 

ro  ro  o 
lo  ro  >«• 


Oj  c  <t  i-H  cr 

o  y  y  I-  3 

or  a  <r  z  I—  CO 

a.  o  o  o  UJ  2  : 

X  u.  u-  o  cr  UJ 

N-  t-  fM 
O  O  fVI 
O  O  fM 


.  a  CD 

:  fM  00 
:  N-'  \ 

:  rsj  O 
u.  o 


I  00  z 
»  2  O 
:  UJ  y 
!  y  y 

HH  o 
I  CO  o 


O  UJ  UJ 

y  00  00 

X  \  V 

<  C  t- 
C»  2  X 
V  Q_  0- 
0-1-1 
>  c  c 


C  nO 

00  t- 


to  3 

O  T-  O 

a  V 

CD  - 

a 


o  o  o 
y  y  y 
y  y  y 
o  o  o 

U)  O 


:  T-  w  h-  ^  ^ 


y  ro 
y  II 

O  l-H 

o  0- 


CO  UJ 
Z  -I 

o 


3  O 
u-  cr 
0- 


p-  Z 

o  o 


O  C  HH 

cr  t— 

.  O  u» 


O  U- 


'  UJ  o 
■  a.  ct: 


X  <r  fr 

<  UJ  i-H 

X  c,  r— 


O  U_ 
UJ  u»  O 


i-H  fr 
00  3 
UJ  •  — J 

r—  09  HH 
c  z  c 

-J  t-H  u- 
3  O 
u)  — I  cr 
— I  UJ  o 


1“  y  2  ♦ 

HH  fr  i-H  ♦ 

<r  o  N^  ♦ 

3  u-  00  «  « 


Q_  O 

cr  2 


UJ  o  u. 

2  Z  3 
•-H  C  00 

r—  r— 

3  00  09 
O  2 
Ql  a  —i 
ro  oo 

3  u_  c 
«  00  O  O  -li 


y  ♦ 
2 


fM  O  ♦ 

\  M 


/-N  o  ♦ 
C  00 

o  o 

I  CD 
Ul  3  ♦ 


O 

vt 


OOUJU)OU)U)UJO 


2 

y 

3 

—I 

o 


63 


Zi 

o 


T-fNJfO'.^m'Oh-oooo  T-f'jfo»'^u>'Or^oo^aT-oj»o»^u^'ON-  oof>»av"fvjKi»^u^'ON-oo^OT-fvjrO'^i/^'Oh-or»^a 

lAlAlAlAirNirNirNU>ir>'O'0'0'0'0^'0'0'0'0f^N-N-N-N-N-N-N-^“N-00  00  00  00  00  00  00  00  Q000^f>.^0'^^0'^O^O 


o  ♦ 


o  * 

f\i 


N- 

o 


I-I  Z  z 

OJ  O  r>i  O 

\  hg  u_  i/>  r>i 

u>  _I  ^  •  _I 

•  O  ^ 

o  O  f-  ^ 

II  O  U-  C  II 

u>  f— •  K'  or  •  _j 
^  ^  a  o 

»—  ro  »—  n  f-  ^ 

<c  ^  w  <£  or  rr  X 

»■  ci  y  <  ru  II 

or  c  or  m 

o  o  uj  o  z  o  y 

u_  Q  Of  U-  X  C3  I— 

c  ^ 

c  o 

o  r: 


^  o 
f-  a 
o  c 

CD  OJ 

\  ^— 

m  II 


O  I 
CD 


a 

I  < 


_j  I—  p^  w  II 
II  i-i  no 
O  W  13 


■»  O  -»  Ll-  > 


rO 
ITS  XV 
•  hO 
CD  XV 

*  • 

♦  O 


«  i-i 

XI  to  tJ 
-»  p^  ♦ 
X  ru  XV  y 
•  y  < 

X  ^  V—  *-l 

D>  CC 

II  n  II  II 

<  y  y  y 

CD  «C  «C  <£ 

TO  i-t  i-t  i-( 

«c  cc  <->  o 


O  /-V  00 

mo 

>■*/-»  a  ^ 

T-  <  •  \ 

'-'CO  O  00 

♦  ♦  I  o  -> 

-rN  ♦  <f  ^  ^  « 

_i  y  y  CD  y  /-v  m 


XV  ^ 
XV  /-N 

o 

XV  oc 
XV  X 
V  * 
O  OJ 


o 
y 

X  ro  I 

>.*  XV  , 


<c  <c  cc 

O'  HH  LLl  < 


I 

y 

>.*  <  < 


«C  X  O  -f 

cn  ♦  y 


f—  <->  ' 

II  II 

a  y  : 


«c  o  o 
X  a  • 
m  +  I 
to  t-  o  : 

OD  I  y  : 

<;  ^  X  I 
'-y  '-y  II  I 
XV  : 


O 

to  m 

*  o 

Q.  O'  /-» 

♦  •  u_ 

OJ  a  to 
'X  OJ  ♦ 
>.*  ♦  /-» 
/-'  <  -J 

y 

_J  X  rvj 


OJ  •  -J 

♦  o  «c 


y  x^  z 
z  y  X 
I  H-  ^ 
X  + 
V-  I  Qi 


/-»  y  y 


I  OJ  ^ 
f  t/)  I 


^<<<c<cu-ii-yo 


O  y-\  /-> 

O  _J  ^  cv  '-y  < 

Xf-t-'x  I  Q'rurja.x 

o  \oytooxv— xz 

rumoJu-H-j'-yX'xLux 

♦  y  II  II  *  II 
oiijoyz<co<c»-‘Q: 

I—  I—  I—  XWCI-V—  li-u-  <c 

HH  II  -I  ^  II  CO 
OaOZ»J-<0<u-Z 
o3tax»-i»j-tD>j-cox 


V-  I 
O  r- 
tO  II 
II  CD 

V-  o 

tJ  O' 

3  a 

r>j 


a 
OU  CD 
OJ  X 


O  -J 


o  'X  'X  c 

V—  A  o  y 
«C  <  QC 


O  UJ 
ta  O' 


Z  z 

o  rJ  '-'  O 

rvj  u.  m  M 

y  /-X  •  y 

^  PO  CD  ^ 

+  II  O  ti-  O  II 

OJ  »-H  X  00  •  _1 

to  %  'X  cr  o 

+  p*^  OJ  V-  II  f-  y 

f-  cr  'X  <£  CK  o  X 

to  P^  o  y  <  ^  II 

M  .  *r  or  X  ^ 

o  UJ  o  y  o  y 


_i 
11 


O'  u.  to  ^  cr  I 


00  p^ 

c  o 


X 


t-  OJ  t- 
X  X  G 
X  X  OJ 


OJ  Oj  OJ 

o 

c 

X 


P^  X 
G  CD 

o-  o 


64 


O 


S  SSSSSSSSSlI^«MMj^rororotOrororororororororororororotOrororoKiroKirorororoKir^KirnK.K.K.K. 


c: 

r^ 


o  ♦ 

>0 


o 

ro 


ro  ro 
K>  ro 

l>0  r<0 

►o  ro 

a  a 

♦  ♦ 


r- 

>0 

>0 

'O 


CD  X 
X 
Q. 


z  a 
^  o 

^  LTV 

a  T- 
o  s 

■>4“  ^ 

%  II 


m  lA 

ru  f\i 
a  o 
c  a 

a  a 

I  I  W  o  ar 

w  w  <  ♦  >- 


* 

I  ^  ^ 
I  ^  < 

Z  Q. 
O  0- 
,  o  «c 
♦  o 


/->  Ci 
_J  o  ^  < 


a.  Cl 

w  XX 

UJ  N.  UJ  UJ 

I-  «-  X  II  11 
l-H  II  ^  ^ 
QC  O  X  O  O 
o  -5  u-  o 


I-  z  ^ 

I—  l-H  II 
UJ  >  < 

II  II  Q. 

z  a.  CL 

l-H  00  « 

>  >  o 


O  -I 
^  Q. 


z  z 

l-H  Z) 
CD  CD 


♦  z 

^  «r 

^  ♦  00 

/-N  X  ^ 

^  -I  z  ^ 

^  w  <  >0 

«c  «c  ^  ^ 

I—  I—  ♦  K 

UJ  UJ  X  >0 

CD  X  Z  — I  ^ 

♦  I-  Z) 

♦  CD 

f-  tA  ♦ 

K  O  X 

>•  o  z 


X 

a. 


I—  ^  w 


> 

•K 

o 

z 

< 

•K 

z 


>  a 


CA 

Z  0. 

O  0. 

O  <  CD  X  «  « 


I-  »-  O  hT 


■n  o  ♦  ^ 
•H  ♦  -K 


^  CD 

w  W  II 

II  M  ^ 
X  X  UJ 

z  o  > 
<  z  a. 

O  <  X 


X  a  a 

I  I  a  NO  ai- 

o.  u.  o  o  UJ  •  «r 

00  u.  O  f\J  Z)  o  ^ 

>  l-H  a  Z  r-  Q. 

II  O  ^  O  IH  '-H  W 

u-  W  II  I—  I—  II  II 

u-  z  Z  Z 

►H  U-  »-  o  O  IH  3 

O  IH  >.  O  CJ  CD  CD 

NO  00 


-J  z 

X  3 

*  CD 

*  * 
CL  Z 
00  l-H 
>  CD 
II  II 
O  UJ 


Kl 
ro 
lA  ro 
•  ro 

O  • 

•K  O  o  c 

♦  -u  < 

^  ^  l-H 

►H  ^  -I  o 


*:  a 

<  * 

»H  « 


o  <  >  »- 


UJ 

^  Q.  X 
<0^10 
I—  u-  <  *>«. 

I  -n  U.  T~ 

Q.  l-H  «  se: 

to  Q£  «4-  >- 

>  >  w  w 

II  II  II  tl 

l-H  o:  s:  s: 

00  <  « 

>  >  UJ  o 


Cl  Q£ 

*  X  T~ 

»:  ^ 

UJ  OJ  CD 
II  II  to 
*!  UJ  II  »-  I-  ’ 
<  I—  I—  Qi  O 
l-H  O  3 
O  >  > 


3 

> 

^  ac 
UJ  to 
I-  > 
O  I 
>  I- 
-I- 

c-  to 
^  > 
II  II 


>o 

a 

f\i 

>0 

o 

f\J 

•  o 

a  ^ 

lA  S 

I  --4  • 

r*N  S  ^ 
^  <4  ^ 
O  ^  ♦ 
U  O  X 

m  o  ^ 
w  a  + 
to  «-  o 
CD  I  s: 
4  X  X 
w  w  II 


to 

♦ 


Cl 

•K 


to  lA 
♦  >0 

l-H  ^ 

CL  O  ^ 
■a  *00 

f\J  C  u. 
f\J 

^  <  _l 

^  s: 


HI  lA  U. 


O  ♦ 

o 


\  -I 
ro  w 


f\j 


a  4 

I  -H  u. 
I  -K  ♦ 
1C  >- 


ro 
ro 

o  ^  ♦ 
y  ro  f\j  z  ar  X 
^  X  ro  I  I  ^ 

•-s.  ro  ^  X  '-H  + 

a  >-v  -I  I-  I  O' 

a  j  ^  0^4 

I  q£ia3clcd 
uja  \Aiartociroi—  X*- 

3f^JU^^Mu.^-HW|OwUJX-5 

z  HI  »■  II  II  ♦  II  to 

i-HOUJOSrZXOXi-HQ£  II 

»— »— r—  I—  x'^“-i— u-  U.HCI— 

Z  l-H  tl  _J  -I  II 


1  lA 

w  ro 
Cl  \ 
X  - 
UJ  T- 

I  - 

It  I— 
X  4 
O  T- 


(O  _J  u.  u.  S' 
>  CD  »-H  IH  >- 


OOQ£0>“ 
OJ  CD  "*  O  X 


u.  <  O  4 


.  O  I 


3  I 

U-  X  o  I 


ro  4  o 

^  3 


'O 

o 

f\i 


4  to  ^ 
ro  ro  a 
ro  4 


w 

3 


s: 

3 


65 


4001  FORMAT(17X,I2,4X,F6.3/4X,F6.3,4X,F6.3,8X,F10.3,17X,F10.3)  351 

222  CONTINUE  352 

RETURN  353 

END  354 


r 

3 


o 

CJ 


66 


<  o 

o  o 

r 

<  LU 

-I  O  • 

o  c 

o  < 

•  Q£ 

O  O 

o 

Q£ 

Q_ 


< 

QC 

O 

O 

oc 

o. 


ro 

oo 

OJ 


o 

OC 


<  o  z 

W  Q  »H 

K 

LU  Lk 

o  o  o 

o 

O  LU 

O  Q£ 

O  C 


<  O 

r  c 


o  o 
a 
o 


z  3r 

r)  o 


O  oJ 


z  o 
o 


o  i 

a: 

a,  to 


u - 1 

O  KH 


h-  c 

C  O' 
>-  Q> 

Cl  ^ 


o 

a  LO 

a 

•  »— 

a  ♦- 

o  o  o 

m  c 


II  ♦- 

o  z 
V  M 

t-  »— 

»—  lO 

Z  M  Z 

LU  ^  l-t  O 
!-♦  Z  U-  l_> 
(_>  LU 


ro 

00  00 
LO  f\i  a 
o  •  O' 
•  a  csi 

a  • 

o 


II 

It  /*>* 

/-%  Z  II 
m  t-i 
z  •  to 
1-1  Z3  z 

•  t_>  r-l 
3  ^  w 

tu  to  to 


•t  -It 

>  -It 


♦-  :r  a  -It 

LU  -K 

o  I  ♦ 

■*c 

u.  to  -Ic 

O  ul  -It 


X  O 
«r  Z 


r3 

o 

to  ♦  rvj 

LU  -K  • 
►H  ♦  N. 


!-•  *-i  X 

"s.  X  to 

*—  -ic 

U.  i_>  O  O 

to  J  LU 

QC  -K 

u-  l-t  »—  <— 

m  ^  z 

LU  « 

LU  u~  t_>  'O 

-j  _i  V 

ft  ■*( 

O  u.  «r  ■ 

W  c  t_> 

O  -K 

tu  LU  u. 

>-  t—  t-l 

cr  * 

o 

1—  u-  X 

Q_  -ic 

to  lU  LU 

i-i  y*  *— 

■ic 

to  _l 

to 

cs  -H 

•f  Z  O 

♦ 

Z  u-  o 

J  * 

y-  o  z 

-K 

LU  o  z 

LU  « 

hH  C  II 

♦ 

O  i-t 

h-i  * 

LU  *—  /-> 

>-  to 

X  -K 

o  o  t—  >- 

♦ 

LU  ^  «C 

to  ♦ 

«£  LU  t_>  CO 

to  >»( 

O  >»«  t_> 

■IC 

X  1-  <  Z 

LU  -It 

X  to  • 

O  « 

LU  O  0,  LU 

t-i  ■¥■ 

<  z  o 

z  « 

^  a:  to 

M 

X  LU  > 

<  -IC 

<3:  0- 

X  ’*( 

t_>  £1  «C 

■n 

LU  -K 

LU  -It 

fV,  ♦ 

J  -It 

»— 

O  ♦ 

r-t  -It 

o:  * 

1—  -It 

n  ^ 

CV.  -K 

tu  -If 

O  1— 

LU  -It 

o  to 

or  “K 

-If 

•  z 

o  + 

o  -n 

o  c 

z  ♦ 

o 

Of  -If 

>0  0  t-> 

o  -n 

o 

CL  -If 

vO  o 

c:^  « 

•  fV- 

-It 

o 

♦ 

*  o  o  o 

OC  -n 

•  o 

-ic 

coo 

O  -It 

o  ^ 

♦ 

oo  •  o 

1-  -It 

II  ^ 

oc  o  • 

fL  -It 

/->  ■ 

+ 

o  ro 

LU  -H 

to  rvj  o 

o 

t->  -H 

z  rvj 

m 

LU  -It 

t-i  00 

II  ro 

Of  -It 

•  II  O  t- 

/-s 

-n 

3  m  O 

to  II 

t->  LO 

CL  /-s  II 

^  •  II 

Li.  to 

to  X  O  os 

w  2  • 

Z  LU  >- 

tU  H-t  to 

►HI-  1- 

>  Z 

<■  LU  I-* 

QC  i-H 

X  y:  > 

»-  O  w 

L9  C  t-i 

Z  ^  QC 

w  Qi  H  )— 

J  LU 

v:  «£  >-• 

0.  U.  to 

to  »—  LU 

O'  Z 

z  LU  x: 

O  >•  O'  LU 

o  y*  c 

»—  ♦—  o  to 

O  l-l 

t_>  r-*  ^ 

Li.  O 

<  >  o 

to  o 

Lc  r-t  «£ 

.  UJ 

♦—  LC  ♦— 

Uj  I-*  Z 

lU  X  •-» 

fL  to  C  < 

Z  L3  to 

-T  Z  LU  »— 

cr 

T  Li,  O'  to 

3  UJ  3 

X 

L0 


<  o 

♦-  a 


X  o 
o 
a 


X 

o 

Z) 

o 

X 

to 

a: 

o 


LU  X 
L3  to 
CK 

C  O 


o 

X 

to 


o 

►H 

to 


< 

L9 

Ui 

z 


<  o 
I-  o 

I- 

QJ  tU 

X  a 
a 
o 


•I  o 

o 

X  c 

LU 

X  o 

CL 

> 

0- 

11 

_1  LU 

«c 

_J  UJ 

/-X 

*t  O 

X 

*t  cr 

to 

a 

o 

z 

o 

LU 

a 

►H 

• 

• 

•If 

o 

00 

-It 

a 

LU 

LU 

«c 

o 

z 

I— 

3 

»-» 

o 

1— 

I— 

LU 

t_> 

LU 

X 

LU 

LU 

1— 

3 

00 

00 

C5 

«c 

z 

OC 

t— 

♦-  3 

►H 

1— 

l-  o 

CL 

Z  O 

z  c 

«c 

OC 

< 

z 

I—  LU 

< 

I—  UJ 

►H 

to  (3 

LU 

to  a 

Z  3 

n 

z  o 

LU 

o  o 

CL 

o  c* 

o 

t_>  • 

< 

t_>  • 

oc 

a 

o 

3 

o 


a 

< 


o 

X 

to 


to 

z 

o 


UJ  X 

o  to 


o 

z 


< 

X 


to  to  z 
o 


67 


CHARGE-T0-"1£TAL  RATIO  Q=  O.SKirfO 

I'JITIAL  VELOCITY  Of  F  R  A  6  viE  NT  S  (  F  P  S)  =  A9>5?.725 

AVERAGE  FRAG'IENT  "lA  SS(  GRAINS)  =  0.1A6E  02 


to 

CT  o  o 

r-  CX3 

o  o  o 

or  •  •  • 

UJ  ^  fO  PO 
►- 
UJ 


d 

ctr 


Z3 

O 


o 

I 


o 


(/>  c  c:  c: 
to  o  o  o 
UJ  rvj  ro 

2  •  •  • 

^  o  o  o 


c 

Ctr 

ol 


30 

o 


K'  N.  r^- 
sT  00 

U_  «—  K'  PO 


f\i  N-  r*- 
>  >C  <5 
u.  OC  O  O 


to 


o  o  o 

tvl  o  K' 
U_  UTS  O 


% 


68 


# 


t 


o  o 

r  I 


j  o 
o 


X 

fj> 

o 


X 

o 

o 

X 

CL 


ro 

oo 

f\J 


z 

o 


z 

o 


:?•  X  uj  -ii 


«c  _t 
z  a 
o 


>■  LA 
»—  • 

-o 

-i 

M  X 
OC  I— 
h-l 

tn  ■« 

o 

X  cc 
Q.  UJ 


Z  UJ 
O  -> 
h-i  o 
»-  a: 
CL 

»— 

3  ST 
CL  ^ 
r-  lA 
O  lA 
UJ  «— 


O 

o 

X 


z 

C9 


O 

O  <A 

O 

•  I— 

O  t- 

o  o  o 

lA  O  S’ 

>»■  rvj 

•  CO 

o 

it  »- 

o  z 

iZ  n  < 

t-  I— 

I-  LC  00 

Z  II  z 
UJ  »—  l-H  O 
M  z  ^  o 
O  UJ 
i—i  i—i  a: 
u.  O  O  O 

U_  l—l  h-  *— 
UJ  U_  O  >0 
o  u_  «t  • 
o  ui  u~  c:^ 

o 

lO  UJ  UJ 

oo 


fO 

oo  00 
lA  rvj  o 
o  •  » 

•  O  'O 

o  • 

o 


II 

II  /-s 
^  Z  II 
oo  l-H  ^ 
Z  *00 
i-»  3  Z 

•  UJ  H-l 
ID  ^  W 

u>  (/>  oo 
'u  CC  oo 
oo  J  UJ 
cc  U^  Z 
^  _1  w 

w  <  UJ 


o 

o 

oo  «  >o 

UJ  «  • 

i-H  *  r^ 
I-  *  ^ 
X  « 

UJ  -n 
CL  * 
o  « 
a:  «  It 
CL  *  ^ 
*  00 
o  -n  z 


o 


•  z 

o  o 

>0  0  UJ 

>0  o 
o 

• 

o 

II  't 

^  t“  • 

to  IM  O 

Z  AJ 

i—i  X 

•  II  O  t- 

3  m  o 

UJ  lA 

^  •  n 

oo  QC  O 

Z  UJ  >- 

►-I  I-  I- 

<C  UJ  H-l 

X  ^  :> 

U3  «r 

£y  II  I— 
^  <r  <r  H-i 
CL  u.  to 

X  z 

o  >■  a:  UJ 

U-  h-  C  to 

UJ  ♦—  I— 

<r  >  UJ  ^ 
u-  i-t  <C 
I-  u_  ►- 
LU  l-H  Z 
Du  00  <  < 
z  UJ  ^ 
X  UJ  q:  </5 
CO  to  <r  z 
O 


o  ♦ 

X  ♦ 
CL  « 
it 

X  * 

o  ♦ 
z  * 
o  •»« 
o  « 

it 


o 
o 
•  >0 

o  rvj  o 
o  ▼-  o 
X  •  c 
X  o  • 

O  'O 

LA 

ro 

li  <4^ 
lA 

(o  n 

(L  II 

u-  to  ^ 
u^  z  • 

UJ  H-t  to 

>  ^  z 

QC  H-t 
»—  to  w 
Z  w  QC 
<  _t  UJ 
I-  <  I- 
OO  I—  UJ 
Z  UJ  z* 

o  z:  <r 

UJ  l-H 

u.  o 

oo  O 


UJ  X  t-( 

Z  to  to 

QC  t-(  ^ 


:  T-  z 

:  o  »-i 


z 

o 

it 

z 

u 

to 

-1 

l-H 

S' 

o 

z 

it 

UJ 

o 

z 

LU 

w 

n 

II 

* 

o 

HH 

IH 

« 

UJ 

UI 

U! 

it 

>- 

to 

X 

it 

_1 

z 

to 

UJ 

I— 

>■ 

it 

LU 

»- 

< 

to 

it 

l-H 

o 

<c 

UI 

UJ 

OO 

oo  * 

to 

l-H 

UJ 

n 

»- 

a: 

< 

z 

UI  it 

X 

oo 

• 

o 

it 

UJ 

UJ 

UJ. 

o 

CL 

UJ 

l-l  it 

<c 

z 

to 

z 

it 

UJ 

-1 

i> 

X 

s 

oo 

1—  -H 

X 

UJ 

> 

< 

it 

—5 

m 

CL 

l-H 

w 

CJC  ♦ 

UJ 

o 

< 

it 

o 

«c 

< 

UJ  it 

UJ 

it 

X 

»- 

CL  ♦ 

it 

CL 

►-I  « 

I-  ♦  z 

UJ  -I*  l-H 
UJ  it 
—>  *  in 
o  *  to 

QC  -H  Q£ 

CL  «  <C 

*  X 
0£  -K  UJ 
O  •»« 

«  u- 
CL  «  O 
UJ  n 
tn  it  h- 
UJ  «  X 
QC  «  to 

♦  HH 

LU 


<  UJ 

to  o 


o 

CO 


«I  o 
»-  o 

I— 

UJ  LU 
m  lA 


o 

X 

oo 


«C  O 
X  o 
CL 

_1  UJ 

«c  o 
X 
X 


z  o 

«c 

»-  UJ 

to  o 
z 

o  >o 

UJ  • 

o 


_J  00 
«t  z 


»-  o 


QC  - 

I-  oo 


> 

«c 

X 


X 

I— 

z 


CL 

CL 


o 

X 

oo 


oo 

z 

to 


z 

<: 


z  o 

c 

«  UJ 
to  X 
rvj 

• 

o 


<  o 
>-  o 
»- 

UJ  UJ 

m  ro 
>0 
o> 


<  o 
X  o 
CL 

_1  UJ 
<r  ^ 
O 


a: 

O 


> 

«c 

X 


o 


z  o 

<  X 

Ul  «c 

to  rvj  ui 

Z  LA  CL 

O  >0  CL 


_J  00 
<  z 


1-0  z 

<:  w  to 

s:  »-  HH 

_l  00 

I-  CL 

UJ  UJ 

to  _l  > 

-  •-» 


oo 


< 

to 


o 

z 


o 

z 


69 


C^ARGE-TO-METAL  RATIO  Q=  0.31160 

I'^ITIAL  VELOCITY  OF  F  R  AGMENTS  (  F  PS)  =  A569.200 

AVERAGE  FRAGi^ENT  M  A  S  S  (  G  R  A  I N  S)  =  0.719E  02 
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CH  ARGE-TO-flETAL  RATIO  3=  0.2S621 

INITIAL  VELOCITY  OF  F R AGNENT S ( F PS ) =  4403.375 

AVERAGE  FRAG^IENT  N  A  S  S  (  G  R  A I N  S)  *  0.3560  02 
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APPENDIX  C 

DIFFERENT  PROJECTILE  CONFIGURATIONS 

Rotational  symmetry  at  the  source  of  explosion  cannot  be 
assumed  when  a  group  of  projectiles  are  in  a  configuration  other 
ttian  the  "square  matrix"  arrangement.  Assuming  that  all  required 
information  is  available,  the  following  steps  would  help  in 
predicting  the  probability  of  detonation  for  grouped  projectiles: 


Region  One 


Region  Two 


■Poiar  angie 


rfK\ 


Region  TWO 


C> 


Region  One 


Arimuthal  angle 


Figures  18  2 

1.  Determine  the  fragment  mass  distribution  in  both 
Regions  1  and  2  (fig.  1). 

2.  Determine  the  fragment  velocities  in  the  same  regions 
also.  These  distributions  are  taken  as  functions  of 
the  azimuthal  angle  only. 


3.  Determine  the  same  distributions  listed  in  1.  and  2. 
above  as  functions  of  the  polar  angle  (fig.  2). 

4.  Resolve  the  distributions  for  each  region  to  obtain 
parameters  FV  and  FM  which  are  now  expressed  as 
functions  of  two  angles: 


FV  =  fi(<j,,  e) 


FM  =  f2((|.,  e) 
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where  <!'  and  6  are  polar  and  azimuthal  angles, 
respectively. 

5.  Determine  C/\,  the  fraction  of  metal  weight  involved  in 
fragment  distribution. 

6.  With  these  modifications,  the  overall  model  can  now  be 
used  to  predict  probability  of  detonation  for  a  group 
of  projectiles  in  any  configuration. 

It  should  be  noted  that  determination  of  the  distributions 
mentioned  in  1.  and  2.  above  can  only  be  achieved  through  field 
tests  . 
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